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Call for submission - Cali for participation 


qCDCq 

Summer Academy 

Recent Research on DERIVE/TI-92-Supported Mathematics Education 

August 25-28, 1999, Gosing, Lower Austria 

This conference is the fifth in a series of conferences on DERIVE/TI-92 -supported mathematics edu- 
cation initiated by ACDCA: Krems (Austria) 1992. Krems 1993, Honolulu (USA) 1995, Saxo 
(Sweden) 1997. 

Typical topics within the scope of the conference are 

• Reports about successful classroom experiments 

• New ways of mathematics teaching and learning 

• Improvements in the design of mathematics software systems based on mathematical experience 

• Examples of successful, new computer-based mathematics texts, lectures, training units, etc. 

• Influence on assessment 

• Similar studies with other computer algebra systems 

How to submit / Time table: 

Please send title and abstract (15-20 lines) to Bernhard Kutzler. 

E-mail address : b. kutzler @eunet. at 

Postal address : Kariglstrasse 5/6, A-4060 Leonding, Austria 

• May 15, 1999: Deadline for submissions 

• May 31, 1999: Notification of acceptance and publication of a detailed program 

• June 25, 1 999: Deadline for early bird registration 

Local and travel information: 

The early bird conference fee is ATS 700 ($ 55) and applies to registrations received by June 25. Af- 
ter that the conference fee is ATS 1000 (S 80). The expected price for a single room with full board 
included is ATS 990 ($ 78,-) per day. 

Gathering point is St, Polten (capital of Lower Austria), from where we will take an ancient nostalgic 
train to the hotel in Gosing. (The train will depart from St. Polten August 25 at 3:25pm.) 

St. Polten can be reached as follows: 

• Plane: Fly to Vienna, then take a train to St. Polten. 

• Train: Take a train to St. Polten from most major cities in the neigbouring countries. 

• Car: Take the A1 motorway towards Vienna, leave it at the St. Polten exit. 

More detailed travel information will be sent out with the detailed program end of June. 

The Alpenkurhotel Gosing, situated at 900 meters sea level 60 km far from St. Polten, is a well hidden 
residence admidst the charming surroundings of the National Park “Otscher Tormauer", eye to eye 
with the huge, impressive Otscher Mountain. Nothing but 15,000 acres of forest around you. The 
hotel offers indoor and outdoor swimming pool, outdoor tennis court, sauna and steam bath, solarium 
and fitness studio. In cooperation with the hotel management we will offer a program for accompany- 
ing guests / family members. 

Conference Co-chairmen: 


Josef Bohm 
Helmut Heugl 
Bernhard Kutzler 


nojo ,boehm@pgv .at 
hheugl@netway . at 
b.kutzler@eunet at 
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Liebe DERIVE- und Tl-Freunde, 

DNL#33 leitet das nun schon 9. DUG-Jahr 
ein, zu dem ich Sie alle recht herzlich beg ru fie. 
An meiner Stelle lacht Ihnen dieses Mai Filou 
entgegen. Sie werden ihm im nachsten DNL mit 
einem Gesprach mit seinem Schopfer JJ Dahan 
wieder und ausfuhrlich begegnen. 


Neben den ewig jungen Themen DGL und 
Fraktale mochte ich Sie besonders auf den Arti- 
kel von Richard Schorn zur Gruppentheorie 
hinweisen. Er bietet uns eine nicht alltagliche, 
aber umso attraktivere Anwendungsmoglichkeit 
eines CAS. Die DERIVE-Implementation kann 
leicht auf den Tl ubertragen werden. 

Es bleibt mir nur noch, Ihnen alien ein Frohes 
Weihnachtsfest und ein gesundes , friedliches 
und erfolgreiches Jahr 1999 zu wunschen. Uns 
wunschen wir, dass Sie wie bisher weiter zur 
DUG halten und uns gewogen bleiben. 

Zum Schluss noch eine sehr wichtige Service- 
leistung von uns gemeinsam mit Soft Warehou- 
se Europe: Sie konnen von nun an alle DERIVE- 
und Tl-files des aktuellen DNLs uber eine neue 
Internetquelle abrufen. Klaus Jurgen Kutzler - 
die nachste Generation - hat auf der Flomepage 
von SWHE eine DUG-Sektion eingerichtet, auf 
der Sie alle files finden: 

www.derive-europe.com/english/dugengl.htm [] 

Damit mussen im DNL nicht mehr alle Listings 
abgedruckt werden. Falls Sie keinen Internetzu- 
gang haben aber gerne etwas haben wollen , 
dann rufen Sie bitte einfach an Oder schreiben 
Sie. Ich schicke Ihnen eine Diskette, wenn Sie 
nicht auf die Jahresdiskette warten wollen. 

Herzlichen Dank Klaus Jurgen fur die angebote- 
ne Hilfe und .... Servus. 

Damit wunsche ich Ihnen vie I Spafi mit diesem 
DNL und im 9. DUG Jahr. 

Bis zum nachsten Mai - und denken Sie bitte an 
Gosingl! 

Josef Bohm, Editor 


Dear DERIVE- and Tl-friends, 

DNL#33 opens a new DUG year. Yes, it is 
number 9. Welcome to you all. This time Filou is 
smiling to you on my place. In the next DNL you will 
meet him again, when he has a interesting TI-92 
discussion with his creator Jean-Jacques Dahan. 

First of all I’ d like to point at the announcement 
of the ACDCA Summer Academy 1999 
on the opposite page and simultaneously 
invite you to come to Lower Austria. 
Unfortunately the Hagenberg Confer- 
ence had to be cancelled but we are sure 
to offer with the Summer Academy in 
one of the most attractive regions of our 
federal state more than only a substitute. 
Because of this announcement I had to 
omit both Book Shelf and web sites as 
well. 

Besides the always 
young topics ODEs and fractals I want 
to point at Richard .Schorn’s group 
theory article. This is a field where - at least in the 
DNL - we do not find so often CAS. The DERIVE 
implementation can easily be transferred onto a TI. 

I am very grateful that Johann Wiesenbauer 
"composed” his TITBITS from a long e-mail discus- 
sion concerning an effective algorithm for calculating 
Stirling numbers. He promised to inform us in one of 
the next DNLs about the use of these numbers. 

Finally I am very happy to announce an impor- 
tant service together with Soft Warehouse Europe. 
You can download all files (DERIVE and TI) men- 
tioned in the actual DNL from SWHE’s web site. K J 
Kutzler - the next generation - has established a DUG 
section on SWHE’s home page where you can find all 
files for downloading: 

www.derive-europe.com/english/dugengl.htm [ * ] 

It is no longer necessary for me to print each 
listing in the DNL - and for you to type it in. If you 
don’t have any access to the Internet then feel free to 
call or to write and I’ll send a diskette, if you don’t 
want to wait for the Diskette of the Year. 

Once more, many thanks Klaus Jurgen for your 
support and cooperation. 

I wish much pleasure with this DNL and for the who- 
le DUG year number nine. 

Hope to meet you in Gosing 

Josef Bohm 
Editor 



Ich mochte Sie gleich auf die 
Ankundigung der ACDCA Sommer- 
akademie 1999 auf der gegenuber- 
liegenden Seite hinweisen und Sie 
herzlichst dazu einladen. Die voran- 
gekundigte Konferenz in Hagenberg, 

00, musste leider abgesagt werden 
Aber wir sind sicher, mit der Som- 
merakademie mehr als nur einen 
Ersatz anbieten zu konnen. Wegen 
dieser Ankundigung muss dieses Mai die Bu- 
cherschau und auch die Sammlung 
der websites entfallen. 



Bon jour, je suis Filou. 


EDITORIAL 
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The DERIVE-NEWSLETTER is the Bulle- 
tin of the DERIVE & TI-92 User Group. It 
is published at least four times a year with 
a contents of 44 pages minimum. The goals 
of the DNL are to enable the exchange of 
experiences made with DERIVE and the 
TI-92/89 as well as to create a group to 
discuss the possibilities of new methodical 
and didactical manners in teaching mathe- 
matics. 

As many of the DERIVE Users are also 
using the TI-92/89 (and TI-NspireCAS) the 
DNL tries to combine the applications of 
these modem technologies. 


Editor: Mag. Josef Bohm 
A-3042 Wiirmla 
D'Lust 1 
Austria 

Phone/FAX: 43-(0)660 31 36 365 
e-mail: nojo.boehm@pgv.at 


Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the DNL. It 
must be said, though, that non-English arti- 
cles will be warmly welcomed nonetheless. 
Your contributions will be edited but not 
assessed. By submitting articles the author 
gives his consent for reprinting it in the 
DNL. The more contributions you will 
send, the more lively and richer in contents 
the DERIVE & TI-92 Newsletter will be. 


Next issue: June 1999 

Deadline 15 May 1999 


Preview: Contributions for the next issues 

3D-Geometry, Reichel, AUT 

Graphic Integration, Linear Programming, Various Projections a.o., Bohm, AUT 

A Utility file for complex dynamic systems, Lechner, AUT 

Examples for Statistics, Roeloffs, NL 

Quaternion Algebra, Sirota, RUS 

Various Training Programs for the Tl 

A critical comment on the "Delayed Assignation" :==, Kiimmel, GER 

Sand Dunes, River Meander and Elastica, The lighter Side Halprin, AUS 

Type checking, Finite continued fractions, Welke, GER 

Share Holder's Considerations using a CAS, Bohm, AUT 
Kaprekar's "Self numbers", Finite Groups, Schorn, GER 
Linear Programming - Graphic Solution on the Tl, Kirmse, GER 
Implicit Multivalue Bivariate Function 3D Plots, Biryukov, RUS 
Physical Problems, Magiera, POL 
LU-Decomposition, Morales, COL 

An Investigation on Calculus with the TO-92, Meagher, AUT 
and 

Setif, FRA; Vermeylen, BEL; Leinbach, USA; Speck, NZL; Biryukow, RUS 
Wiesenbauer, AUT; Aue, GER; Koller, AUT and others 


Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wiirmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag. Josef Bohm 
Herstellung: Selbstverlag 
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Jean-Jacques Dahan, Toulouse, FRA 

Dear Friend, 

I supposed you are interested by publishing into DERIVE letter the article I’ve promised to you at the 
end of the last year. I send the translated (in English) version in order that a great number of readers 
can read it. I am working on a Spanish version. 

DNL : Many thanks for your charming, enchanting and inspiring contribution about your talks 
with your little friend Filou. As you can see he took my place in the Letter of the Editor and he 
said that he wouldn't leave that position until he hasn't seen his conversation with his master 
published. 

Julio Cesar Morales, Universidad Nacionai des Colombia, COL 

Dear Mr. Boehm, 

I send you the LUACTR.MTH file along with the corresponding .DOC and .DMO-files to see if they 
qualify for inclusion in your newsletter. LU ACTR.MTH includes the functions U_P(A) and L_P(A) 
which (unlike the function that once appeared on the DERIVE web site) give the LU factorization (LU 
decomposition) of matrix A even if row interchanges are necessary. 

DNL: Dear Mr. Morales, 

As Editor I am feeling very happy and proud to have - after a long time of absence - an active 
South American member in the DUG and your article will be printed in our bulletin in the next 
future. Muchas gracias. 

Gunter Scheu, Pforzheim, GER 

Dein Artikel fiber die "Logos, " hat mir ausgezeichnet gefallen und in Deinem Artikel "Das Pas- 

cafsche " hat mich sehr gefreut und sogar fiberrascht was der TI-92 doch alles kann. Ich habe fiber 

160 Seiten fur mein Buch geschrieben und hoffe, dass es in diesem Halbjahr noch beim Dfimmler- 
Verlag erscheint. 

Der Unterricht mit dem TI-92 macht mir und meinen Schfilem enorm Freude. Fast in jeder Stunde 
leme ich noch was Neues. 

Now in English: I have two questions concerning the TI-92 (Standard Model): 

© Why is the TI-92 unable to calculate the integral [ (x(x - l)(x - 2 ) 2 dx ? 

© There are also problems in calculating the 1 st derivative of f(x) = ^ f x - ^\dx for x = t. 


Calculate the zeros of the first derivative and then substitute them into / without the ABS. Why 
does the TI present the solutions as ABS-expressions? 
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DNL: Dear Gunter, many thanks for you friendly words about the last DNL. I really appreciate 
the comments of experienced teachers and "advocates of the GAS" - to quote Carl Leinbach. 
To your questions: As you can see the integral works, but it needs an incredible amount of 
time to find the integral value. Changing the upper bound to 1.9 or 2.1 accelerates the proc- 
ess? 

I sent both requests to Michelle Miller - my very reliable link to the Tl-experts - and via her I 
received a very short answer from the experts: 

1. and 2. are TI-92 bugs that were fixed on the TI-92+. 

There was a third answer from Michelle, because I added one more problem. We - Wolfgang 
Propper and I came across this bug (???) when we worked on our book about integration 
methods. In one section we dealt with the representation of integral functions. This is a prob- 
lem on the PLUS, which doesn’t occur on the Standard Model. Compare two TI-92 screen 
shots given below: 

— X JOT X < 1 

<=> when(x<1 ,-x,x) -> f(x) 

x for x>\ 

® All integrals work - with one exception: a general upper bound: 


Standard TI-92 TI-92 PLUS 




f FE* Y F3¥T 'i 

R I gebralcal c [Other 

f FS ’| 

PrgfilO 

FS<r V 

Cl tan Or! 

■ f ** ^ 1 * R(>A 


Done 

r-x,x< i . 
lx, else 

SCf txi-x-JLt > 

MAIN 

RAD EXACT 

FUNC 2/3© 


The Tl's answer: 3. occurs on the 89/92+ because of a change that was made in symbolic 
integration due to a problem report. 

My conclusion: This is another example for the fact that there is no advantage without any 
accompanying disadvantage. 


See how TI-NspireCAS performs: 






D 
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It is interesting to see how DERIVE6.10 - and other CAS - are coping with this integral: 


# 1 : 


fCx) := 

If x < 1 

-X 

X 


Using the IF-construct gives a wrong result. 
DERIVE simply substitutes the boarders into the 
both functions! 



# 5 : gCx) := - x-xt-™, x, 1 ) + x-xCL x, «0 

5 27 

# 6 : f g(x) dx = 

-2 2 


# 7 : 


t 

I g(x) dx 
0 



'l 


1 1 
-SIGNCt - 1 ) 


2 


2 


Taking the CHI-function works correct. See be- 
low the integral with two symbolic boarders. The 
TIs refuse to give any result. 


b 

r 

1 

2 ^ 
a 


r 2 > 

b 1 


\ g(x) dx = 
a 

< 2 

2 ^ 

-SIGNCa - 1 ) + 

.2 2 j 

-SIGNCb - 1 ) 


I recommend reading the very interesting article contributed by our Canadian friends 
Genevieve, Michel and Henri Integration of Piecewise Continuous Functions in DNL#91. 



► CAS 0 

1 

Q 

f(x):=lf[X<1,-X,X] 

f ( x ) : = If [1 > x, —x, x] 

2 

m 

\ i ? i ?1 

Integral: If 1 > x, — - x\ - x £ + ci 

3 

O 

lntegral[f(x), x, -2, 5] 

27 

-+■ — 

2 

4 

O 

lntegral[f(x), x, 0, t] 

-+ ? 
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wxMaxima & CASIO ClassPadll 


(%il) 

f (X) 

:=if x < 1 then 

-x 

else x; 

(*ol) 

f(x): 

= if x< 1 then 

- X 

else x 

(«i3) 

wxplot2d ( [ f (x) ] , [x, 

-5, 

5] ) $ 


(%t3) 



(%i4) integrate (f (x ) , x, -2, 5); 

5 


(%o4) | if x< 1 then -x else x dx 

- 2 


define f(x) = 


-x, x< 1 
x, x>l 


f (x)dx 


-2 

□ 

f (x)dx 

□ 


done 


13. 5 


/I 


-x,x< 1 

X, x£l 


dx 


Mathematica 


In [14]:= f[x ]:= Piecewise [{{-x, x < 1}, {x, x^ 1}}] 


In [16]:= f [X] 


Out [16]= 


-x X < 1 
X X > 1 


In [17]:= Integrate [f [x] , {x, -2, 5}] 

27 

Out [17]= — 

In [18]:= Integrate [f [x] , x] 

x<\ 


Out [18]= 


2 


-1 True 


In [19]:= Integrate [f [x] , {x, 0, t}] 

-x x < 1 


Out [19]: 


Out [26]: 


'■r 

= Int 

'=r 


x> 1 


dlx 


In [20] := Integrate [f [x] , {x, y, to}] 

, — x x < 1 


. titx 


x> 1 

In [21]:= Plot[f [X], {x, -5, 5}] 



Peer van de Sanden, NED 

Hello, I want to compute Stirling numbers of the first and second kind: 

The following recursive functions work but are slow ... 

DA/L: This email (you will find it on page 23) was the beginning of a very interesting discus- 
sion on generating Stirling Numbers. Johann Wiesenbauer dedicated his TITBITS column to 
this topic (page 23). 


Riideger Baumann, Celle, GER 

I send a little “Knobelei” (= Brain Teaser) - puzzle for the DNL-readers which shows DERIVE ' s ca- 
pabilities to work with lists: 

What is the rule for the sequence: 1, 11, 21, 1112, 3112, 122213, 312213, 212223, ...? 

Rudeger’s original file from 1999 does not work with DERIVE 6.10. I had to replace his orog- 
inal definition of H(v,k) (see below) by a SELECT -construct which is expression #1 in Rude- 
ger’s hint file. Josef 

DIMENSION^) 

H(v, k) := I IFtv = k, k) 

1=1 i 
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#1: H(v, 10 := DIM(SELECT(a = k, a, v)) 

Rtideger’s Hint: 

#2: COUNTS (v) := APPEND(VECTOR(IF(H(v, k) > 0, [H(v, k), k ] , []), k, 1, 9)) 

#3: COUNTS ([1, 1, 2, 3, 3, 4, 1, 2, 4, 3, 3]) 

#4: [3, 1, 2, 2, 4, 3, 2, 4] 

#5: COUNTSEQCvO , n) := ITERATES (COUNTS (v) , v, vO, n) 

DIMENSION(v) - n 

#6: POT(v) := VECTOR (10 ,1,1, DIMENSION(v)) 

#7: ZF(vO, n) := VECTOR(w-POT(w), w, OOUNTSEQ(vO, n)) 

#8: New simplify: 

#9: ZF([1], 10) 

I believe that it would not be so easy to program this sequence e.g. in PASCAL. On the other hand the 
Hofstadter Sequence 

1, 3, 7, 12, 18, 26, 35, 45, 56, 69, ... 

can easily be produced by a PASCAL program. Who of the DERIVE family can help with a DERIVE 
program? 

DNL: Dear Riideger, I can imagine that your request is a challenge for some of the DNL- 
readers. (Among others I count on Johann Wiesenbauer.) By the way, I know your Teaser- 
Sequence in a slightly different variation: 

1 , 11 , 21 , 1211 , 111221 , 312211 , 13112221 , 1113213211 , ... 

In the very recommendable CRC Concise Encyclopedia of Mathematics ( Eric W. Weisstein) I 
found a nice and for me unknown name for this sequence: “Look and Say Sequence”. 


The Hofstadter Figure-Figure Sequence is defined as follows (also from CRC): 

Define F(1 ) = 1 and S(1 ) = 2 and write: 

F(n) = F(n-1) + S(n-1), 

where the sequence {S(n)} consists of those integers not already contained in {F(n)}. For 
example, F( 2) = F(1 ) + S(1 ) = 3, so the next term of S{n) = S( 2) = 4, giving F( 3) = F( 2) + S(2) 
= 7, so S(3) = 5 and F(4) = F( 3) + S(3) =12. Continuing in this manner gives the “figure” se- 
quence as shown in Rudeger’s problem. 

(Hofstadter, Godel, Escher, Bach, p. 73 ( Seite 79 deutsche Ausgabe) 

Sebastiano Cappuccio, Forli, ITA 

Dear Derivers, 

VECTOR (IF(x<0,2x+k) , k , - 5 , 5 ) works well with DERIVE 2. 08 and I can simplify and plot it, 
getting a sheaf of parallel half-lines. 

Why does it neither work with DERIVE 3.01 nor with DfW 4.112 How can I get the same sheaf with 
them? 

Thank you, Sebastiano 


p 
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Peer van de Sanden, NED and Carlos Fleitas, Comenar Viejo, ESP 


Hello, I don't know why it does not work in DfW 4.11 but I think that it has something to do with 
nesting formulas. This will work: 

F(x,k) : =IF (x<0, 2x+k) 

VECTOR ( F (x, k) , k, -5 (5) 



John Alexiou, USA 

The problem seems to be that simplification does not occur for the 2 nd or 3 rd arguments of the 
IF( ) function for example , an expression 


#1 : IF (u>d, u-d, u+d) 
when applied to 

#2: LIM(#l,d,t) 


you would expect the result IF (u>t,u-t,u+t ) but the simplified result is 

#3: IF(u>t,LIM(u-d,d, t) , LIM (u+d, d, t ) ) 

no matter how many times you simplify it does not simplify the 2nd and 3rd argument of IF( ). 

I always use the STEP( ) or CHI( ) functions which work better, unless you call STEP(O) which re- 
turns 1/2 +/- 1/2 which cannot be evaluated. 


The problem exists for: vector ( if () ) , lim(if()), and iterate (if () ) and possibly more ... 


Soft Warehouse, Inc. (Albert Rich) 

When DERIVE simplifies an IF statement, it simplifies the "conditional” clause (i.e. the first argu- 
ment) of the IF statement. If its truth- value can be determined, the "then" clause (i.e. the second argu- 
ment) OR the "else" clause (i.e. the third argument) of the IF statement is simplified, as appropriate. 
Note that the other clause is NOT simplified. Therefore, if the truth-value of the "conditional" clause 
cannot be determined, neither the "then" clause nor the "else" clause is, or should be, simplified. 

For example, if x is unknown and IF(x<0, 1+1, 2+2) is simplified, DERIVE does not simplify 1+1 or 
2+2, since the truth- value of x<0 is unknown. 

However, the arguments of user-defined functions, as distinct from the clauses of IF statements, ARE 
always evaluated before calling the function. Therefore, you can define replacement functions for the 
IF statement that do evaluate the "then" and "else" clauses as follows: 

IF_THEN (condition, then) : =IF (condition, then) 

IF_THEN_ELSE (condition, then, else) : =IF (condition, then, else) 
Then, to generate the vector of IF statements that you want, you can simplify the command 
VECTOR (I F_THEN (x<0, 2x+k) , k, -5, 5) 

Hope this helps clarify the situation. 

Aloha, Albert D. Rich 


D 


-N-L#33 


DERI VE & 77- 92- USER - FORUM 


P 9 


There was another discussion concerning "Isolated points", raised by an experience made 
by one of James Eckerman's students. 

James Eckerman. Sacramento, California 

A student of mine asked why her TI-89 gave her a graph that has points in quadrant two for the func- 
tion cos ( sqrt ( x ) ) when she had the mode set on "REAL". I got the same picture with DERIVE. 
I'm sure some- one out there knows the correct explanation. Please help us out! Thanks. 



{ rin-V F2-" F3 " FH T F5-V FU- "FF-# iTl 
]T$$l5|Zem Tk4« Re€raph]M<ith)Di'aw Pm 'Z j 





MAIN Mtf EKfl 1 

:t func 


Harald Lang, SWE 

Look at the Taylor expansion of cos(.v): 

6 4 

X X 

TAYLOR (GOS(x), x, 0, 6) = + 

720 24 



If we substitute sqrt(x) for x, then we get: 

6 5 4 3 2 

X X X X X X 

TAYLOR COOS G/x) , x, 0, 6) = + + +1 

479001600 3623300 40320 720 24 2 

where the right hand side makes perfect sense also for negative values for x. I suppose what you see is 
the right hand side above. 

Adri van der Meer, NED 

e x _|_ e ~ x 

Because cos(/x) = cos(-/x) = , and this is real if x is real. 


If you choose the ’’real” mode in DERIVE , you don't prevent it from using complex calculus for "in- 
termediate results". You only force DERIVE to choose the real branch (if possible) in the case of mul- 
tivalued functions. 
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Peer van de Sanden, NED 

Hi Jim, here is a simple function that also gives a plot for negative x-values: 

F (x) : = RE (x + SQRT (x) ) 

In this case it is expected. 1 am asking for the real part of a complex number. The output is real al- 
though the part in the brackets (x + SQRT (x) ) is complex. 

Your function COS (SQRT (x) ) does the same thing. 

The real mode works on the input and on some simplifying algorithms. It doesn't affect the inbetween 
results. In your case: 

SQRT() COS() 

x(real) —> z(complex) — » y(real) 

Look at the following DERIVE simplification for the complex numbers a + b i and b i and compare 
the results: 

-b + L-a b - L - a. b -b 

e e e e 

[G05(a + b-O, C0SCb-O] = + , + 

2 2 2 2 

The second expression has no imaginary part in its simplification! The output is real. 

r e^-x e -^~ x 

Therefore (for x < 0) : cos( v x ) = — — H — - — . 

The following function gives the expected plot: COS( i f ( X>0, V( x) ) . 

Johannes Wiesenbauer, Vienna, AUT 

James example reminded me of another one which is still a mystery to me: 

Why doesn’t DfW 4.1 1 show any points in the third quadrant when plotting the curve SQRT(x) A 2? 
Anyone out there who has a clue? (It goes without saying that "dirty tricks" like simplifying this ex- 
pression are strictly forbidden!). 

Cheers, Johann 

J. H. Frisbee, Houston, USA 

The SQRT(x) A 2 problem may have to do with the way roots and powers are computed. First, it is in- 
teresting to note that if one uses the product of SQRT(x) * SQRT(x) the curve is complete. Now sup- 
pose SQRT(x) is replaced by EXP(LN(x)/2). The more original form EXP(LN(x)/2) A 2 exhibits the 
incomplete character and, unlike the above product form EXP(LN(x)/2) * EXP(LN(x)/2) is also in- 
complete. 

I have no idea how DERIVE actually handles the LN and EXP functions but I suspect that they-are 
probably used for powers (and roots) and that is the more fundamental cause of the incomplete plot. 


James Eckerman, Sacramento, California 

When I read James Eckerman's email, I remembered a strange function, which I obtained last fall. 
I tested my just bought TI-89 trying to calculate the locus of the maximum curvature of the power 
function curves. 
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In the course of that calculation I plotted the function / (x) 


x — 2 

x 2 (2x-l), 


and was rather aston- 


ished to see isolated points within a gap. Further work with DERIVE led me to some interesting re- 
sults as you can see with the enclosed MTH-flle: LOCMAXC.MTH (which is among the files in 
MTH33.zip). 




Compare the plots: TI-89 left above, Voyage 200 others. You can see the “isolated points”, 
Lorenz is writing about, Josef. 

As I mentioned in the MTH-file-comments the reason for the occurrence of the isolated points is the 
negative sign of the basis combined with integer values of the exponent of the power term. The sim- 
plest example for this phenomenon is the (mostly complex valued) function F(x):=(-l) A x or 
"squeezed” in F(x):=(-l) A (l/x). Compare the well known F(x):= x A x, They all produce spirals in the 
complex space (x, RE(y), IM(y)). 

This is the spiral in an 
isometric representation 
(using ISOMETRIC from 
GRAPHICS. MTH - DOS- 
Times, Josef). 


Ray Girvan, UK 

I haven’t looked through your derivation in detail; but I did try checking the plot of/(*) against 
GrafEq , which is normally very good at ‘difficult’ functions with singularities. GrafEq suggests only 
one singularity (1.5, -1/9) in the gap between x = 0.5 and x = 2. 
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Johann Wiesenbauer, Vienna, AUT 


Hi all, the function 

F (x) := ( (x-2 ) / (x A 2 (2x-l) ) ) A (1/ (2 (x-1) ) ) 
introduced by Lorenz Jaeneke is a very remarkable one, indeed, though some things are far less mys- 
terious when viewed from the right angle. (Hence, those who love mysteries are advised to skip the 
rest of his note!) 


When drawing its parametric representation in the complex plane, namely [RE (F (x) , IM (F (x) ] (by 
the way pros will input this term in the form [ ref x, imf x] and leave it to DERIVE to decipher it!) 
for x in [0.5, 0.99] and [0.99, 2] (adjust a suitable range before plotting, e.g. [-0.12,0.12] both forx 
andy) you’ll see at once what is going on here: If x approaches 1 then its image f(x) will move on 
one of two spirals that both have the same circle around the origin as limiting curve. Its radius is 
given by 

LIM(ABS (F (x) ) ,x, 1) = e A (-5/2 ) 

Those "isolated points within a gap” are simply the countable many intersections of the two spirals 
with the real axis. In particular, they are certainly not singularities, as claimed by Ray Girvan (though 
I like his homepage very much!), because there is only one genuine singularity of f(x) when viewed in 
the complex plane, namely x = 1. 1 hope you all will be enchanted by the picture of the two spirals as I 
was and there will be no troubles with older versions of Derive. 


Needless to say that I f m always referring to the current version of DfW which is 4.11 
at the moment. If you have an older version DfW 4.xx instead use the option of a 
free update on Derived homepage! Johann 

Cheers, Johann 


Ray Girvan, UK 

Very neat. Another useful view is to plot RE ( F ( x ) ) and IM ( F ( x ) ) separately, which shows F(x) 
as out-of-phase real and imaginary sinusoidal components, both increasing in frequency as they ap- 
proach x — 1 . The "isolated points" are where IM = 0 and RE has maxima and minima. (In this case, it 
looks as if GrafEq was on the wrong track). 

One thing I notice: when I plot RE (F (x) ) I get a smooth plot. But if I apply Simplify/Basic to 
RE(F(x)) and plot the result, the section of the curve between 0.5 and 2 becomes discontinuous. Any 
idea whaf s going on. 

I am sorry! I know: singularities are places where the function blows up to infinity, not isolated solu- 
tion points. I wasn't thinking. Thanks for the clarification. 


^ ^ ^ ^ aMM A 2 ^ x 


X 

t AN, /I A l\ ^ 

0.5 2 


v v fJfv X x 

0.5 1 1.5 Z 


Plot of RE(F(x)) 


Plot of IM(F(x)), {Josef 1999) 
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Plotting real and imaginary parts now: 



Soft Warehouse, Inc (Albert Rich) 

Hello Ray Girvan 

The result of simplifying RE(F(x)) produces an expression mathematically equivalent to RE(F(x)) but 
expanded out. When the expanded expression is evaluated, roundoff error produces a small imaginary 
component in the result that prevents it from being plotted. 

The moral of the story: If you wrap RE around a complex 
expression in order to plot the real part, do not simplify it. 

Ray Girvan, UK 

It also looks quite good with an offset to show it from a different viewpoint: 

[x+RE (F (x) ) , IM (F (x) -x) ] , (setting the parameters to [-3,3] plots the function reasonably in 
one operation). 

I like the way an apparently discontinuous function like F(x) makes far more sense in complex space: 
just a smooth curve that twirls into imaginary over parts of its range. 


Johann Wiesenbauer, Vienna, AUT 

Hi all, 

I'm very sorry, but the second interval for x should be of course [1.01,2], otherwise you won't see the 
second spiral... As I just tried out it will also work, if you draw both spirals at once by choosing the 
interval [0.5,2] for x, though I prefer the first option because it shows the limiting circle between the 
two spirals more distinctly. 

By the way, isn't it strange that Derive draws this highly complicated curve without batting an eyelid, 
but has troubles with the simple function sqrt(x) A 2? Well, actually Derive is on the horns of a di- 
lemma: In the latter case simplifying the expression before plotting it would fix the problem, but in the 
ease of our spirals simplifying the parametric expression above is not exactly advisable as Ray has 
noticed ... 

On this occasion, let me point out that Forenz did a great job, when computing the actual values of the 
isolated points, namely 


p 
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1 h ( m a nonzero integer) 

2m 

so we'll look over the enormous troubles he had when computing the two real limit points of f(x). Just 
one remark that might be helpful to throw light upon matters of this kind: When computing the limit 
of an expression like 

r \ 

8 /7? 2 ( 1 - Am) 

((2/w + l)(4m + 1) 2 ) 2m 

you have to bear in mind that m is viewed as a complex variable by Derive (no matter what adjust- 
ments had been made before!), so Derive was perfectly right in saying that this limit doesn't make 
sense. It took me some time to find out why Lorenz "succeeded" in finding a "limit" all the same: It's 
because In x and ln(abs(x)) have the same derivative, so when applying L'Hospital's rule the outcome 
is the same as if you had taken the limit of the absolute value of the expression above which does ex- 
ist! A nice fallacy, don't you think so? 

Cheers, Johann 


lim 

m— >oo“ 


You can see the enchanting picture of 
the two spirals - Johannes 
Wiesenbauer - together with Ray Gir- 
van’s curve “twisting into imaginary”, 
Josef. 

I cannot reproduce the plots with 
DERIVE 6.10, Josef. 



Two more “generator - generated” “Koch like curves (page 50) 
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A Spigot - Algorithm for n and e 

Peter Witthinrich, Liibeck, Germany 

In 1995 Rabinowitz and Stanley [1] published a new algorithm for computing the decimal digits of the 
irrational numbers n and e. It uses no floating point routines but only integer arithmetic and produces 
the digits one at a time without using them after computing. So there was no need to store them in an 
array or string [5]. The authors called it 'spigot algorithm' (called "Tropfelverfahren" in [3] - in an 
exercise and example for Delphi programming). ( For those of you who have never heard the word 
spigot - in German it is a “Zapfhahn ”, Prost! Josef) 

The authors explicitly refer to a method published by Sale [2] in 1968 - nearly thirty years earlier, for 
computing the digits of e, the basis of the natural logarithm. It is very short, and it also pumps the 
digits of e one at time, again without using them after computing. Therefore the use of arrays for stor- 
ing the already computed digits, i.e. in [4], can be avoided. 

The algorithm is based upon the well known series of e 



The sum right of 2 is <1. 1 shift the decimal point, i.e. I divide by 10 and simultaneously I multiply by 
10 each of the single summands. I obtain the following decimal place from the integer part of the re- 
mainder within the parentheses. Again we find a fraction as remainder: 



The digits of the intermediate results are calculated by integer division and stored. The last carry de- 
livers the next decimal digit. Algorithmically you need a loop: 

Q“>carry 

For j , $ , 2 , -1 

coeff [ j ]*10+ca rry^temp 
intDi vtteinp, j ) scarry 
rema 1 n ( temp t j )->coeff [ j ] 

EndFor 

Output row, col umn T 1 Part ( carry ) 



The algorithm in [1] was written down in Algol 60, a good structured language, and it is not very dif- 
ficult to translate it into TI-92 programming language (filename: spigote.92p): 
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The output on the TI-92 display has been formatted to 270 decimal places, nine rows with 30 digits 
each. Computing time for the 270 decimal digits was about 35 minutes. The computed digits for e are 
correct compared to the result in [4], Unfortunately the corresponding algorithm for ji[ 1] is more com- 
plicated and longer than the algorithm for e. Both programs are rewritings of the programs given in the 
references in the TI92 language, they both don’t include new own ideas. 

A result of 270 decimal places here needs about seven hours computing time, the TI-92 is too slow: 


spi gote( ) 

Prgm 

Local digits.s, test, row, col umn , i , j , coef f , car ry , temp 
Cl rIO 

Request "Humber of decimals: ", digits 

exprtdigits )-*digits 

4-»s 

( di gits+ 1 )* 2 . 30258509 ^test 


While s+(ln{s)-l)+.5*ln(6. 2831 B52+s latest 
s+l-»s 
EndWhile 
i nt C s )-»s 

Output 0, 0, "Euler 1 s number: e = 2." 

setMode! "Exact/Approx" ."Exact") 

10-s-row 
0-*col umn 

seq(l,n,l,s,l )-*coeff 
For i ,1 .digits 
0-»carry 
For j , s , 2, -1 

coef f [ j ]*10+ca rry-=»lemp 
intDiv(temp, j)-*carry 
remain! temp, j )->coeff [ j] 

EndFor 

Output row, column, iPart(carry) 
col umn+8-»col umn 
If col umn>239 Then 
O^col umn 
row+10^row 
End If 




TTrr Y rs Y ~ 

£ I £ |S dPraniOlSE -i .5=” :: , 


Eulersch# Zahl * & a 2* 

7 1828 1 828459 045235 36028747 1352 
662497757 247B936 9 99595 7 496696? 
627724876630353547584571382178 
523166427427466391932603059921 
8 174 13596629 043572 900334295260 
595630738 132328627 843480763 233 

8 298 S3 7 53 1 95 25 1 i 19 01157363418? 
930702154889149934884167509244 
761460668032264800 1684 774 1 1 353 

M HI H DEGftETFiDii PAR 1/J4 



^ V :: v y f^FiT V 



Kreiszahl ff = 3. 

141592853589793238462643333279 

5 0 28841 9 716 9 3993751 0 5 8' 2 0974944 
592307816406286208998628034825 
34211 70679821 48Q86513282306647 
093844 60955058223 1 7 2 5 359488123. 
481 1 17450284102701 938521185559 
644622948954930381 964428018975 

6 6 5 9 3 3 4 4 6 1 2 S 4 7 5 6 4 8 2 3 3 7 8 € 7 8 3 i 6 5 

2.71^l^g-lJ-^_o.3 £6 9,2 3 , 4 6 0 3 4 9 6 10 
NifiM 6E6»?Pftmi Pft¥ 


End For 

set Mode ( " Exact /Approx” /'Approximate”) 
EndPrgm 


spigotpiQ is included in MTH33.zip. 
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Subject: Shading with DERIVE 

Peer van de Sanden, Netherlands 

Dear Derivers, 

I started this new year (happy new year to you all) with creating an utility file that I needed often in 
the past but never had the time to make it. It's called shade .mth. This file contains four functions 
for shading an area between the two functions/^) and g(x) on the interval [< a , b]: 

SHADE (a, b, space) 

SHADE_UNDER (a,b, space) 

SHADE_BACK ( a , b , space ) 

SHADE_BACK_UNDER (a,b, space) 

You must define the functions X*) and g(x) before using the functions. The parameters a and b are the 
boundaries of the interval [a, b] (with a < b). The parameter space is the space between two spacing 
lines. This parameter is optional. The standard value for space = 0.2. The vertical boundaries x = a and 
y = b will not be calculated and therefore will not be plotted. 

The SHADE -function draws lines in the area where fix) > g(x) on the interval [a, b]. The angle of the 
lines is 45° if the x-scale and the y-scale in the plot are equal. 

The SHADE_UNDER-function draws lines in the area where g(x) >/(x) on the interval [a, b\. The angle 
of the lines is 45° if the x-scale and the y-scale in the plot are equal. 

The SHADE _BACK-function draws lines in the area where f(x) > g(x) on the interval [a,b\. The angle 
of the lines is -45° if the x-scale and the y-scale in the plot are equal. 

The SHADE_BACK_UNDER-function draws lines in the area where g(x) > f{x) on the interval [a, b]. 
The angle of the lines is -45° if the x-scale and the y-scale in the plot are equal. 

The shading works with various functions, The functions do not work on intervals on which either /(x) 
or g(x) has an asymptotic behaviour. 



As I produced a shading tool some times ago (DNL#1 5) and G P Speck presented his " Con- 
centric Curve Shading" in DNLM31, you may think that there is some demand on tools to 
beautify plots. (W.Propper's contribution in this DNL works in the same intention). I had some 
ideas to improve Peer's useful tool and contacted him how to add the boundaries , to com- 
bine both shadings and how to produce a SHADE_BETWEEN-function. After a very short 
while he wrote back: 
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I tried several times to develop an easy to use shading-tool. The functions you added make it possible 
to make life even more easier. Here are my suggestions: 


• Include the boundaries in my original user functions. Then they can be deleted from your func- 
tions. 

• Rename the functions in a more logical way 

• Add one more which shades normally and backwards between the two functions /(x) and g(x) 


Then we will have the following functions: 

User functions (not frequently used): 

SHA_A (a, b, space) 

SHA_U (a, b, space) 

SHA_AB (a, b, space) 

SHA_UB (a, b, space) 

SHA_ABO (a, b, space) 

SHA_UBO (a, b, space) 


User functions (frequently used): 

SHADE (a, b, space) 
SHADE_BACK ( a , b , space ) 
SHADE_BOTH (a,b, space) 


Most people will only need your functions. Thus these functions should have the most comfortable 
names. The SH_A- and SH_U-function look a little like the SHADE-option on the TI-92 

Note: in the Options Menu on the plot screen you must select Points-connected and Points-small and 
perhaps no color cycling. 

Note: I remember that I had problems with functions with asymptotic behaviour like /(x) = l/x A 2 with 
g(x) = 0. Derive will not be able to simplify function SHADE (-1,1). With space = 0.2 the maximum 
will be at x = 0 and you will get a vector with an infinite amount of elements. In this case an other 
value for space will do e.g. SHADE ( - 1 , 1 , 0 . 3 ) . 


I wanted to have shading functions for plotting histograms (Dutch word, I don't know if the English 
name will be the same). I have included an example in histogram. mth. The original shading- 
functions are adjusted for this application. 

You are welcome to publish it in the International DERIVE & TI-92 Newsletter. You additions de- 
serve that the original MTH-file will be adjusted, preferable in the way that I suggested. 

But, would it not be nice if Derive had a shading option or an area color option in the plot screen like 


the TI-92 is providing? I would like to right click 
ground color or shading for that area. 

Groeten van Peer van de Sanden 



k] := [0, 100, 200, 250, 300, 400, 500] 
freq := [20, 40, 50, 75, 50, 10] 
PLOT_HISTOGRAM(k 1 , freq) 


the mouse on an area and choose my own back- 



F(x) :i 0.3-Cx + 2. 5) - (x + l)-(x - 2.7) + 1 
2 

G(x) := 0.05-0 - 2-x) • (x + 4-x + 5) 
SHADE_BOTH(-3.5, 2.8, 0.15) 
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See some two examples: 




f 

f 71 ^ 

3 ^ 

AreaBetweenCu rves 

SINOO, 2 -COS 

2 - x 

, x, - 2 -tt, 71, y 



L 4 J 

2 J 



The next plot was created using the DERIVE 6 option to transfer the plot for a short 
while to the very comfortable Graphics Utility Paint. Click on the plot after copying it in 
the Algebra window and open the Edit Menu for proceeding. 


v , Unbenannt- Paint 


□ 


Start Ansicht 

% Ausschneiden 
"l Kopieren 


| 1 


Einfugen 

Zwischenablage 


tz£ Zuschneiden 

£3 GroBe andern 
Sk Drehen * 

Bild 


! /<>■ ; 

Tools 


||\^ODOZiA'- 
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Umriss - 
& Fullen 



Strichstarke 





Farben 
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Graphic - Tools for the TI-92 

von W. Propper, Niirnberg 


1. Einleitung 

Nach dem Entpacken der Datei GRA_TOOD.EXE haben Sie auBer dieser Dokumentation 
RA_TOOD.DOC eine TI-92 Gruppierung GRATOOD.92G erhalten. Mit Hilfe von TI-Graph-Link 
kann dieses group file in 5 TI-92 PRGMs (*.92P files) entpackt werden. Diese konnen dann (auch 
unter Verwendung von Graph-Link) in den TI-92 iibertragen werden. (Es empfiehlt sich, vorher mit 
2nd [VAR-LINK] ein eigenes Verzeichnis auf dem TI-92 anzulegen - Z.B. UTY - die Programme 
dorthin zu laden und in diesem Verzeichnis ablaufen zu lassen.) 

Die Programme wurden von W. Propper erstellt und werden zur nichtkommerziellen Nutzung fur den 
allgemeinen Gebrauch freigegeben. Obwohl die Programme mit Sorgfalt angefertigt wurden, kann 
eine Gewahrleistung, gleich welcher Art, nicht ubernommen werden. 

Verbesserungsvorschlage, Anregungen und Kritik werden dankbar (und moglicherweise zahneflet- 
schend) entgegengenommen. 

Snail-Mail: W. Propper, Josef-Simon-Str. 59, 90473 Niirnberg 

E-Mail: w. proepper@ wpro . franken . de ; 

2. Programmbeschreibung 

Der TI-92 erlaubt schnell und einfach Graphen von Funktionen darzustellen. Will man diesen Graphen 
jedoch ein vorzeigbares AuBeres geben, stoBt man schnell an Grenzen. Die Plot Styles ’’Dot”, ’’Square” 
und ’’Thick” stellen, wenn iiberhaupt wirksam, nicht das Optimum an Sauberkeit dar und die iiber die 
Funktionstaste F7 im Graphik Fenster erreichbaren Linien- und Textwerkzeuge dienen in den wenigs- 
ten Fallen dazu, einen Graph mit aussagekraftigen Zusatzen zu versehen. Die fiinf Programme der 
GRA _ TOOLs Gruppe sollen hier eine gewisse Abhilfe schaffen. Mit einem der Programme konnen 
Graphen mit verschieden geformten Markierungen fur Punkte versehen werden und mit den vier rest- 
lichen Programmen lassen sich strichlierte Linien zeichnen. 

(Short explication in English follows at the end of the contribution, Josef) 

pntat(x,y,p,l) zeichnet Punkte mit den Koordinaten (x,y). Dabei konnen x undy sehr unterschiedliche 
Formen im Aufruf annehmen (Beispiele s. unten). Die Form der Markierung wird mit dem Para- 
meter p bestimmt. 1st p = 1, so wird ein ”+”-Kreuz und bei p = 2 ein ”x”-Kreuz gezeichnet. Einen 
kleinen Kreis erhalt man mit p = 3. Bei p = 4 entsteht ein gefulltes Quadrat (■) und bei p = 5 ei- 
ne gefiillte Raute (♦). Der Parameter I bestimmt die GroBe der Objekte. Bei I = 2 betragt ihre ma- 
ximale Ausdehnung 5 Pixel, bei I = 3 ist sie 7 Pixel. Nach regularer Ausfuhrung des Programms 
schaltet der TI-92 in das Graphik Fenster. Liegen jedoch alle Punkte auBerhalb des Plotbereichs 
oder wenn fur p bzw. I unzulassige Werte angegeben wurden, bleibt der Home Screen aktiv und 
der Programmaufruf meldet done (obwohl eigentlich nichts passiert ist). 

Beispiele: pnt(1 ,2,1,2) erzeugt ein kleines ”+”-Kreuz im Punkt (1,2), wahrend pnt(2, 1 ,2,3) bei 
(2,1) ein groBes ”x”-Kreuz zeichnet. pntat({-1 ,2,3},{1 ,-1 ,2},1 ,3) setzt groBe ”+”-Kreuze an die 
Stellen (-1,1), (2,-1) und (3,2). 
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Wenn die Funktion/(x) deklariert ist und bei x = 2 eine stetig behebbare Definitionslucke besitzt, 
dann setzt pntat(2 Jimit(f(x) 5 x,2) 5 3,2) an diese Stelle einen Kreis mit Radius 2 Pixel. 

Um alle Punkte mit horizontalen Tangenten (in den meisten Fallen die Extrema) mit einem gro- 
Ben M + M -Kreuz zu markieren, ruft man pntat(zeros(d(f(x) 5 x) 5 x) 5 f(zeros(cy(f(x) 5 x) 5 x)) 5 1 ,3) - oder 
besser: zeros(c/(f(x),x),x), was eine Liste der Nullstellen der 1. Ableitung berechnet, und dann 
pntat(ans(1 ),f(ans(1 )),1 ,3) - auf. 

Das Markieren aller Nullstellen einer Funktion mit einer groBen Raute wiirde mit 
pntat(zeros(f(x),x) 5 0,5 5 3) allerdings nicht funktionieren, da der x-Parameter (bei mehr als einer 
Nullstelle) eine Liste mit mehr als einem Element liefert, der Wert 0 fur den ^-Parameter jedoch 
nur als einelementige Liste interpretiert wird. 

Mit f(zeros(f(x),x)) anstelle von 0 erhalt man jedoch das gewiinschte Resultat. 

slv(x,l,p) zeichnet eine vertikale strichlierte Linie an der Stelle x. Die Lange der Striche ist I, die 
Zwischenraume haben die Lange p. Da ein Funktionsgraph durchaus mehrere vertikale Asympto- 
ten besitzen kann, darf der Parameter x auch eine Liste sein. 

slh(y,l,p) zeichnet entsprechend eine horizontale Linie bei y. Die Parameter I und p haben die glei- 

che Wirkung wie bei der vertikalen strichlierten Linie. (Fur y ist nur eine Zahl bzw. ein numerischer 

Ausdruck zulassig.) 

sls(x1,y1,x2,y2,l,p) erzeugt eine strichlierte Linie durch die Punkte (x7,y7) und (x2,y2) mit Strich- 
lange I und Zwischenraum p. Diese Prozedur wird normalerweise fur schiefe Asymptoten oder 
auch Wendetangenten verwendet. 

sis ruft slv bzw. slh, wenn die x- bzw. die y-Koordinaten gleich sind aus dem Folder uty. Falls die 
Programme in einem anderen Folder des TI-92 installiert sind, fiihrt dies zu einem Programmfeh- 
ler. 

slf(Up) liefert schlieBlich eine strichlierte Naherungskurve der Funktion y = t{x). Dazu muss t ein 
Ausdruck mit der unabhangigen Variablen x sein. Die Naherungskurve besteht aus Geraden- 
stiicken, deren Anfangs- und Endpunkt auf der Kurve y = t liegen. Die angegebenen Parameter fur 
Lange und Zwischenraume beziehen sich auf horizontale Stiicke und sind bei einem Koordinaten- 
system mit gleichen Langen der x- und der y-Einheiten (z.B. ZoomDec) stimmig. Bei starker 
Verzerrung (z.B. ZoomStd) sind die fur den Betrachter sichtbaren Strecken gleich lang, obwohl, 
bezogen auf das Koordinatensystem, die steileren Stiicke langer sind als die flacheren. 

Der Term t darf auch langere Liicken im Defmitionsbereich haben (z.B. V(x 2 - 1), der in ] -1,1 [ 
nicht defmiert ist). Allerdings lauft slf dann einigermaBen langsam. Um „passende“ Werte fur I 
und p zu erhalten muss man ggf. etwas experimentieren. 


Die fiinf Programme laufen in den 4 zweidimensionalen Graph Modes. Es empfiehlt sich, zuerst den 
Graph der Funktion zeichnen zu lassen und anschlieBend die oben gezeigten „Versch6nerungsarbei- 
ten“ vorzunehmen. Das Kommando ReGraph (F4) im Graph Screen loscht alle von den oben genann- 
ten Funktionen erzeugten Linien und Punkte. 
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Ein Polynom 3. Grades: Extrempunkte mit 
Kreuzen, Wendepunkt mit "x"-Kreuz, Nullstellen 
mit Rauten. Wendetangente strichliert. 

Cubic: Turning points inflection point "x", 
zeros tangent in inflection point dashed. 



Eine gebrochen rationale Funktion mit einer verti- 
kalen Asymptote und einer Parabel als asymptoti- 
scher Kurve. Besondere Punkte sind wie vorher 
gekennzeichnet (Nullstelle mit "■") 

Rational function with a vertical asymptote 
and a parabola as asymptotic curve. Zero is 
marked by a 



Eine gebrochen rationale Funktion mit vertikalen 
und schiefen Asymptoten. Punkte sind gekenn- 
zeichnet wie vorher. 

Die stetig behebbare Definitionslucke ist mit ei- 
nem Kreis markiert. 

Rational function with two vertical and one 
slant (oblique) asymptote, points are marked 
like on the left graph, the removable discon- 
tinuity is marked by a circle. 



Eine Kurvenschar (f k (x) = (-2 k + In x)-ln x), in 
der die Kurve der Tiefpunkte kurz gestrichelt, 
die der Wendpunkte mit langeren Strichen ein- 
gezeichnet ist. 

Family of curves (f k = (-2k + In x)-ln x), with 
locus of turning points dashed (short seg- 
ments) and locus of inflection points (long 
segm.). 


Summary of the tools 

pntat(x,y,p,l) point at (x,y), p for the mark (1 ,2, 3, 4, 5) for + , x,#, ■, ♦ ; 

I for size: 2: small, 3: large 

slv(x, I, p) vertical dashed line at position x, with segments of length I and spaces of 
length p 

slh(x, I, p) horizontal dashed line at position x, with segments of length land spaces of 
length p 

sls(x1, yl, x2, y2, 1, p) segment between points (x1,y1) and (x2,y2), length and space 
slf(t, I, p) dashed function graph for function t(x), length and space 
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Titbits from Algebra and Number Theory(14) 

by Johann Wiesenbauer, Vienna 

Back again! As some of the readers may know there has been a vivid discussion about Stirling 
numbers of the first and second kind via the DERIVE-mailbase quite recently. In the aftermath 
Josef has asked me to sum up the results of this discussion and to tell you more about Stirling 
numbers in general, in particular how they are defined and what they are good for. Well, actu- 
ally this stuff belongs to combinatorics (in fact, their DERIVE-implementations can be found 
in the new utility- file COMBINAT.MTH !), it is true, but I have no guilty conscience at all 
when treating this topic here, as Stirling numbers also have a number of important applications 
in algebra and number theory. 

The discussion mentioned above was triggered off by the following problem posed by Peer van 
de Sanden on the DERIVE-mailbase(derive-news@mailbase. ac.uk): 

Hello , 

I want to compute Stirling numbers of the first and the second kind. The following recursive func- 
tions work but are slow: 

Stirling numbers of the first kind: 

st-irllCn, k) := 

If n = k 
1 

If k = 1 
Cn - 1)! 

Cn - 1) - sti r"ll(n - 1, k) + stirllCn - 1, k - 1) 

Stirling numbers of the second kind: 

sti r!2(n , k) := 

If n = k 
1 

If k = 1 
k 

k -sti rl2(n - 1, k) + stirl2(n - 1, k - 1) 

For this last function I have found a faster one: 

n p 

k COMBtk, p) - (k - p) -C-l) 

stirl_2(n, k) := I 

p=0 k ! 

Does anyone know a function for the Stirling numbers of the first kind which is faster then the one 
from above: 

Here are the recursive rules: 

Stirlingl(nJ) = (n-1)! 

Stirling 1 (n,n) = 1 

Stirlingl(n,k):= (n-1) * Stir ling 1 (n-1, k) + Stir ling 1 (n-1, k-1) 

Peer van de Sanden 

I named the function stirl because later we will use STIRLING (name of the respective 
DERIVE functions). Johann will give a short description what the St. n. are good for. In this 
revised version HI take the space to add some more accurate notes on them, Josef 
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and comparing the coefficient of x k for 0 < k < n , the recursion formula for Stirling numbers of the first 
kind mentioned by Peer is easily obtained. 

My answer contained the following passage: 

.... To the best of my knowledge in the current version of the new utility file combinat.mth the fol- 
lowing formulas are used for Stirling numbers: 

lk k - j_ n 

stirling2(n, k) := I COMB(k, j_) • (-1) -j_ 

k! j_=l 

stirlingl(n, k) := 

If n < k 
0 

If n = k 
1 

(-l) A (n + k) •I((-l) A j_-C0MB(n - 1 + j_, n - k + j_) • 

C0MB(2-n - k, n - k - j_) • sti rl i ng2(n - k + j_, j_) , j_, 1, n - k) 

The implementation of STlRLlNGl(n,k) above is certainly faster than that of Peer van de Sanden 
(by a factor 2 or so). Since it essentially follows a suggestion of mine , I would like to pass on the 
question: Is there still a faster one? ... 

Yes, it is true, A1 Rich had asked me exactly the same question so me months before and I had sug- 
gested in my answer the implementation of STIRLINGl(n,k) above (however, without considering the 
cases n < k and n = k\ by the way, Peer too overlooked the case n < k and also the case k = 0 in his 
implementation). Since my implementation was so much faster on average than Peer's (speaking of a 
factor 2 only was actually a colossal understatement!), I thought that my final asking whether there 
was still a faster one was only a rhetorical question. Little did I know then, how much I was mistaken! 

This is where the third protagonist, Ralph Freese of Hawaii, comes into play with the following mes- 
sage: 

Here's a faster Stirling 1 version : 

sti rl i ngl_aux(v, n) := VECTOR(IF(k = 1, n!, IF(k = n + 1, 1, n-v + v )), 

k k - 1 

k, 1, n + 1) 
n - k 

stirlingl(n, k) := (-1) • (ITERATE(sti rl i ngl_aux(v, DIM(v)), v, [1], n)) 

k 


Notes: 

1. This is only valid if 1 <= k <= n; it really should have an IF statement like Derive's version in 

front of it. 

2. I've included the sign as Derive did; Peer van de Sanden 's numbers are called unsigned Stir- 

ling numbers of the first kind. 

Concerning calculation time: 

Peer van de Sanden algorithm is slow because to calculate Stirling 1 (2n,n) two calls are made to 
stirlingl(2n-2,n-l), four calls are made to stirlingl (2n-4,n-2), ... , and 2 n calls are made to Stir- 
ling 1(1,1). So this implementation requires exponential time. 
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The program above is really the same recursion as Peer van de Sanden's except that it remembers 
what it has previously calculated. 

The Derive program is much better: it has running time proportional to n"3. The program above 
has running time proportional to n/: 2. 

Shortly after he added the following Post Scriptum. 

The final n should be n - 1 and I noticed that Derive was returning the unsigned Stirling 1 so the (- 
l) n+k should be removed. Stirling 1 (n,k) should be defined: 

sti rl i ngl(n , k) : =ITERATE(sti rl i ngl_aux(v , DIM(v)) , v, [1], n - 1)) 


Well, should I really have to swallow the bitter pill and accept that I had clearly failed in finding the 
best way of tackling this problem? At first, I wasn’t prepared to do this as you may conclude from my 
following reaction: 

Hm, ... although Ralph ’s implementation of Stirling numbers of the first kind is certainly streets ahead 
of anything I have seen so far on that score, it remains to be seen whether it is really faster than mine 
in practice (though it should be asymptotically faster, if the estimates above are correct). 

Unfortunately, when trying to check this I came across unexpected difficulties. For small n, say 
n<100, my implementation is notably faster, it is true, but for larger n (n=103 and k = 1 seems to be 
the smallest example) the function SUM( ) in stirlingl (n,k) 

lk k - j_ n 

stirling2(n, k) := I C0MB(k, j_) • (— 1) -j_ 

k! j_=l 


stirlingl(n, k) := 

If n < k 
0 

If n = k 
1 

(-l) A (n + k) •I((-l) A j_-COMB(n - 1 + j_, n - k + j_) • 

C0MB(2-n - k, n - k - j_) • sti rl i ng2(n - k + j_, j_) , j_, 1, n - k) 

all of a sudden goes to pieces returning senseless values such as fractions. 

To see that my routine should work and SUM () is actually the culprit all you have to do is to replace 
SUM( ) by SUM(VECTOR( ...)) in stirlingl (n,k), but it goes without saying that this remedy is sheer 
madness when it comes to performance. 

Thus in my opinion in terms of speed a fair comparison between the two implementations seems to be 
impossible before a fix of this bug and has to be postponed ...On the other hand, it may well be that 
very few Derivers, if any, are troubled by this ... 

(By the way, according to Albert Rich this SUM error was simply caused by a dummy variable name 
conflict: stirlingl and stirling2 both used the same dummy variable j_ ...) 
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This is what Ralph answered: 

Absolutely right. (Referring to the first paragraph of my email!) 

I didn ’t get the SUM error. On my computer computing Sl(2n,n) was faster with your version than 
mine even with n=100. But the time as a function of n did seem to be growing at a slower rate than 
yours. Viewing Sl(n,k) as a lower triangular matrix, may method calculated the first n rows and re- 
turned the (n,k) th entry. But it really only needs the entries in the parallelogram bounded by (1,1), 
(k,k), (n,k), (n-k+1,1). 

Here is a program to do that: 


stirlingl_next2(k, row, shift, offset, u) := [row, offset + shift, 


VECTOR(IF(j = 1, (row - l)-u , IF(j = row, 1, (row - l)-u 

1 j - offset + 1 

+ u )), j, offset + shift, IF(row < k, row, k))"| 

j - offset J 

sti rl ingl_nextl(v, k, n) := stirlingl_next2(k, 1 + v , IF(v > n - k, 1, 

1 1 


0), v , v ) 
2 3 


stirlingl(n, k) (ITERATE(stirlingl_nextl(v, k, n), v, [1, 1, [1]], n - 


D) 

3,1 


sti rl ingl(10, 5) 


269325 


(Sorry, this isn ’t very readable.) The area of the parallelogram is k-(n-k+l). The worst case is k = n/2. 
it is only about twice as fast for k=n/2 but for k near 1 or n, it is much faster. Yours is worst when k is 
small since it sums from 1 to n-kso this one is faster than yours when k is small. 

Ralph 

When checking Ralph’s arguments above it dawned on me at last that he was right and I had to change 
my basic strategy. This is from my answer: 

(...) After looking into the whole business once more, I must admit that you have a point there, on top 
of leading to that SUM-bug the formula I used is also very inappropriate when it comes to small val- 
ues of k. Therefore I have made up my mind to change sides and try to beat you at your own game. 
What about the following implementation? 

Sl(n , k):=IF(k>n , 0 , IF(k=l, (n-1) ! , IF(k>n/2 , (ITERATE([DELETE(n_- APPEND(v_, [0])+ 

APPEND([0] , v_)) , n_+l] , [v_, n_] , ITERATE ( [n_* APPEND (v_, [0])+APPEND([0] , v_) , 
n_+l] , [v_,n_] , [[1] ,0] ,n-k) ,k))aU, (ITERATE([DELETE(n_. APPEND([0] ,v_) + 
APPEND(v_, [0])) ,n_+l] , [v_, n_] , ITERATE([n_- APPEND([0] , v_)+APPEND(v_, [0]) , 
n_ + 1], [v_, n_] , [[1], 0], k), n - k))aU))) 
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S1C10, 5) = 


269325 


S1C100, 2) = 


4831850496458679681600550566304058591440571607234875011918788469532628073581063287428- 

459115887885635393993680616672360632116716255156633600000000000000000000 

S1C1000, 999) = 499500 

And here are some execution times on a Pentium 116 PC (1999) 

SI (1000, 2): 0.5s vs. 2.3s 

SI (1000, 5 00): 20.2s vs. 161.9s 
SI (1000, 999): 0.5s vs. 2.2s 

Now (in 2014) are the execution times: 0.063s (2566 digits number), 2.27s (1627 digits num- 
ber) and 0.063s for the 6 digits number. Josef 

I hope this top performance will reconcile you even if you don ’t agree with my programming style due 
to the missing auxiliary functions. (Sorry!) 

Cheers, Johann 

I’d like to explain the decisive idea behind my implementation above, because it could turn out useful 
in your own programs. Similar to Pascal’s triangle one can compute the “ Triangle of Stirling Numbers 
of the first kind ”, e.g. 

VECT0R([n, VECTOR(Sl(n, k), k, 0, n)], n, 0, 3) 

'0 [ 1 ] 

1 [ 0 , 1 ] 

2 [ 0 , 1 , 1 ] 

3 [0, 2, 3, 1] 

4 [0, 6, 11, 6, 1] 

5 [0, 24, 50, 35, 10, 1] 

6 [0, 120, 274, 225, 35, 15, 1] 

7 [0, 720, 1764, 1624, 735, 175, 21, 1] 

.3 [0, 5040, 13063, 13132, 6769, 1960, 322, 28, 1] . 


Now it is obvious that it should be much faster to compute each row of this triangle as a whole vector 
by using the built-in vector operations than computing it element wise. Take the row number 5 for 
example: Instead of computing each of its components 

0 = 4-0 + 0, 24 = 4-6 + 0, 50 = 4-1 1 + 6, etc. 

we can compute it at once from two modified copies of row number 4 (zeros have to be inserted in its 
beginning and end, respectively) in the following way: 

4-[0, 6,11,6,1,0] + [0, 0, 6, 11,6, 1], 

This accounts for the frequent occurrence of expressions such as 

n_-APPEND(v_, [0]) + APPEND([0], v_) 
in my implementation above. 
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Only later did I exploit Ralph’s idea that only the Stirling numbers of the first kind in a certain paral- 
lelogram are actually needed to the full extent. It is contained in my following email to all Derivers 
connected to the mailbase (skipping many other interesting contributions in-between in order to cut a 
long story short and come to an end): 

Just in case you also belong to the people who dot the i \s and cross the t ’s like me, you might enjoy the 
following a little bit tidier ’version ’ of my STIRLING l(n,k): 

sti rlingl(n,k) := 

IF(k > n, 

MAX(1 - k + n, 0), 

IF(k, 

0, 

IF(k = 1, 

(n - 1)!, 

IF(k = n - 1, 

C0MB(n, 2), 

IF(k > n/2 , 

(ITERATE([n_-DELETE(v_) + DELETE(v_, 1 - n_ + n) , n_ + 1], 

[v_, n_] , ITERATE([n_.APPEND(DELETE(v_) , [0]) + v_, n_ + 1] , 

[v_, n], ITERATE ([n_- APPEND (v_, [0]) + INSERT(0, v_) , n_ + 1] , 

[v_, n_] , [[1], 0], n - k), 2-k - n) , n - k)) , 
(ITERATE([n_-DELETE(v_, 1 - n_ + n) + DELETE(v_) , n_ + 1], 

[v_, n_] , ITERATE ([n_-v_ + APPEND(DELETE(v_) , [0]), n_ + 1] , 

[v_, n_] , ITERATEf [ n_ • INSERT (0 , v_) + APPEND(v_, [0]), n_ + 1] , 
[v_, n_] , [[1], 0], k), n - 2-k), k))j,lj,l))))) 

Numerical experiments seem to indicate that the asymptotic growth of computation times for Stirling 
numbers of the first kind is much better than that for Stirling numbers of the second kind, though the 
latter can be given explicitly by the simple formula 

lk k - j_ n 

stirling2(n, k) := I C0MB(k, j_) • (-1) -j_ 

k! j_=l 

This leads to the following interesting question: Is this inevitable or does exist an implementation for 
Stirling numbers of the second kind as well that is asymptotic faster than the one above? 

Shortly after I gave the answer myself: 

Yes, there is an implementation with a far better behaviour and as a matter of fact it is so simple that 
it is almost embarrassing. (Why couldn ’t this guy do his job properly from the very start?!) 

sti rling2(n, k) := (ITERATE([APPEND(v_, [0]) + INSERT(0, v_) , n_ + 1] , [v_, n_] , 
[[1], 0], k))|l-VECTOR((-l) A (k - j_)-j_ A n, j_, 0, k)/k! 

For small n and k, you may gain some fractions of a second by using the old implementation, but in 
my opinion it doesn ’t pay off to split the computation depending on the size of the arguments ... 

Again it was Ralph who found a flaw (though harmless this time!) in the routine above thereby putting 
the finishing touches to it: 

Nice! 

I did note that you didn ’t need to have the n_ since it is not used. So it can be simplified a little to: 

sti rling2(n, k) := ITERATE (APPEND(v_, [0]) + INSERT(0, v_) , v_, [1], k) • 
VECT0R((-l) A (k - j_)-j_ A n, j_, 0, k)/k! 

Ralph 
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Well, there may still be minor improvements, e.g. 

sti rl i ng2(n , k) := IF(n < k, 0, ITERATE(APPEND(v_, [0]) + INSERT(0, v_) , v_, [1], 
k) • VECT0R((-l) A (k - j_)-j_ A n, j_, 0, k)/k!) 

but its performance is already staggering in the present form. For example, try to compute sti r- 
1 ing2 (2000 ,1000) with DERIVE and other CAS (nomina sunt odiosa!) and you’ll see the differ- 
ence, believe me!. The same goes for sti rl i ngl above, where I recommend comparing the computa- 
tion times for sti rlingl(1000,500). 

Note: (sti rling2(2000, 1000) is a 3355 digit number, which needs 1.41 sec calculation 
time, sti rling2(1000,55) gives a 1627 digit number needing 0.93 seconds.) 

Concluding another nice programming task: sti rlingl(n,k) can also be defined as the number of 
permutations on n letters having k cycles, whereas sti rl i ng2(n , k) is the number of ways to parti- 
tion a set of n elements into k disjoint subsets. How could DERIVE-routines for sti rl i ngl(n , k) and 
sti rl i ng2 (n , k) look like that make use of this? 


Additional comments on Stirling numbers: 

I found many interesting explanations, proofs and identities in one of the “computer science 
bibles” - Concrete Mathematics by Graham, Knuth & Patashnik. 

Here is a small selection together with the DERIVE verification: 


[ 1 ] 

[ 2 ] 

[3] 


n 

-H 

k 

14 

I-' 

\=—i 




m\ m\f^\k ) 
m + n + 1 

m 


integer k, n 
k"(- \) m ~ k 
Y,(n + k) 


k=m 


n + k 
k 



M 

n 

[n + ll 

k 

[4] 


=2> 

1*4 




k = 0 

m 


(-ir 


[1] does not work with the functions given in Titbits14, but it works with the DERIVE - func- 
tions. 


[STIRLING1C10, 7), STIRLING2(-7, -10)] = [9450, 9450] 

stirling2(13, 7) = 5715424 

17 13 7 - k 

I C0MB(7, k) -k -(-1) = 5715424 

7! k=l 


19 

stirlingl(19 +8+1, 19) = I (8 + k) •stirlingl(8 + k, k) 

k=l 


60383004803151030 = 60383004803151030 

20 11 - k 

I stirling2(21, k + 1) -stirlinglCk, ll)-(-l) = 167960 

k=ll 


C0MB(20, 11) = 167960 
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A finite Group with linear rational Functions 

Richard Schom, Kaufbeuren, GER 
richard. schom@t-online.de 

Expressions A and B (see below) can be used to constmct a finite group. The binary operation 
P<x)Q PQ is defined by replacing variable .v in expression Q by the expression P which is “P in Q”. 

A=— , B = 3-x 

x — 1 

Find the order of the group by building the complete group table. Explain your work using DERIVE. 


Solution: 

CaseMode := Sensitive 
x 

A := , B := 3 - x 

x - 1 

B 

3-x 
3 - A 

1 

2 

x - 1 

2-x - 3 

x - 1 


Let’s find AB 

Replace x in B by A (A in B) 
and factorize! 

this is now AB ! 


0P(E1 , E2) := lim E2 
x->El 


Introduce function 0P(E1 , E2) for “El in E2” 


OPI(el, e2) := ITERATE(e2 , x, el, 1) 


1 

OP(A, B) = 2 

x - 1 


same for purists who don’t like the limit 
at this occasion. 


1 

OPI(A, B) = 2 

x - 1 


2-x - 3 

FACTOR(OP(A, B)) = once more AB. 

x - 1 


OP(A, A) = x 


The powers of A: A 2 

and B 2 (does not need DERIVE). 


OP(B , B) = x 
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FACT0R(0P(el, grp)) = 


2-x - 3 x - 3 2-x - 3 


x, 


, 3 - x, 


x - 1 


x - 1 x - 2 


x - 2 


FACT0R(0P(e2 , grp)) 
FACT0R(0P(e3 , grp)) 
FACT0R(0P(e4, grp)) 
FACT0R(0P(e5 , grp)) 



2-x - 3 x - 3 ' 
x - 2 ’ x - 2 . 

x 2-x - 3 ' 

x - 1 x - 1 . 

x - 3 

, x, 3 - x 

x - 2 

2-x - 3 x 

x - 1 ’ x - 1 . 


FACT0R(0P(e6 , grp)) 


2-x - 3 2-x - 3 


x - 2 x - 1 



X 


This is the result: 


, = X, A = - 


-,B = 3-x, AB = ^, BA = iA 


ABA- 


x - 3 
x - 2 


with among others BAB = ABA, ABAB = BA, BABA = AB, . . . 


See the operation table: 



E 

A 

B 

AB 

BA 

ABA 

E 

E 

A 

B 

AB 

BA 

ABA 

A 

A 

E 

AB 

B 

ABA 

BA 

B 

B 

BA 

E 

ABA 

A 

AB 

AB 

AB 

ABA 

A 

BA 

E 

B 

BA 

BA 

B 

ABA 

E 

AB 

A 

ABA 

ABA 

AB 

BA 

A 

B 

E 


On my request Richard Schom sent three additional examples: 

o A = — which gives a cyclic group of order 5 (see next page). 

2(x + 1) 

o A=— , B=l-x gives a group which is isomorph to the group constmcted by 


x 


2x-3 
x - 2 


3 - x 


(= group of cover mappings of an equilateral triangle). 


2 

2-x 


o A = 2 - x, B = 


leads to a group of order 8. 
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I tried to create the group table belonging to Richard Schorn’s first additional example. Cy- 
clic groups are groups with all elements being powers of one single element. They can con- 
sist of finite or of infinite many elements. 

Using the TI-NspireCAS I found out the five elements of this group: 

A _ ^ + 3 2 _ (V5+3)(jc + 1) 3 _ (>/5+3)(2jc->/5-1) 

’ 2(x + l)’ " 2x -a/5-1 ’ " 2((a/5 + 1)-x + 2-(a/5 + 2))’ 

a 4 = (a/5-1)-((a/5+1)-x + 2-(a/5 + 2)) _ -2x-V5-3 

4x 2x 


op(e l t e2):= factor/ lim (e2)\ ► Done 


\x el J 


-{J 5 + 3 ) -(Vi" + 3 ) 


2-x+2 2' (x+l) 


( \ -(^+3)- (x+l) 

aa:=op(ajaJ ► 

2-x-JS -1 

( ^ (j5+3){2x-j5-l) 

aaa-oplaa.aj - ( , 

2- ((j5 +l)-.x+2- (,/5 +2j) 


, , -«r-i) ((jr + i) x +2 (^ 2 )) 

aaaa:=op(aaajaJ ► 

4' x 

f v "(Vi - — 2 )- (2' y?+2)'X+5' Jb+ll) 

op(a,aaaJ ► 

2' x 


opfa.aaaj-opfaaa.a) ► o/K 
op(aaaaja) ► x 
op(a.aaaa) ► x 

n 

" 

0 


And this is the operation table of this commutative cyclic group: 



U 

A 

A 2 

A 3 

A 4 

u 

U 

A 

A 2 

A 3 

A 4 

A 

A 

A 2 

A 3 

A 4 

U 

A 2 

A 2 

A 3 

A 4 

U 

A 

A 3 

A 3 

A 4 

U 

A 

A 2 

A 4 

A 4 

U 

A 

A 2 

A 3 
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ON THE SOLUTION OF A LINEAR DIFFERENTIAL EQUATION OF 
THE ORDER n WITH CONSTANT COEFFICIENTS 

Emilio Defez Candel , Vicente Soler Basauri 

Departamento de Matematica Aplicada 
Universidad Politecnica de Valencia, P. O. BOX 22.012. Valencia. Spain 

1. Introduction: 

In this paper we study the solution of the homogeneous linear differential equation of nth order with 
constant coefficients, not covered in Derive utility files, by implementing the file: 

CONST.MTH 

which enables the general solution of the homogeneous linear differential equation of the order n with 
constant coefficients: 


a n y n) ( x ) + a n-i y n ~ l) ( x ) +•••+ ay(x) + a 0 y(x) = 0 . 

The functions in this file are also used to solve the linear differential equation of the nth order 

a n y n) ( x ) + a n-\ y n ~ l) ( x ) +•••+ a x y'(x) + a Q y(x) = f(x) , 
by means of the undetermined coefficients and the constant variation methods. 

The organization of this paper is as follows. In section 2 we review the theoretical results on linear 
differential equation of n th order (see Flett). In section 3 we study the implementation of the file side 
CONST.MTH used for the solution of homogeneous linear differential equation, illustrating its use 
with several examples. In section 4 we apply the CONST.MTH file to the solution of the complete 
differential equation, depending on whether we use the undetermined coefficient or the parameter 
variation method. The file was implemented with Derive XM. 

2. Linear differential equation of n th order 

Definition 2.1 A linear differential equation of nth order is a differential equation that can be written 
like this: 

a n (x)/ ) (x) + a„_i( x )/ _1) (x) +...+ a,(x)/(x) + a 0 (x)j(x) = /(x). (1) 

We can write the equation (1) in a more compact way: 

Ay = f{x) , A = Y j a i (x)D‘ , (2) 

i = 0 

d d n 

where D is the derivative operator, defined as D = — , D n y = — — , D°y = y . 

dx dx 11 

The existence and uniqueness of local solutions of equation (1) is guaranteed by the following result: 
Theorem 2.1 Let all the functions a t (x) , i = 0,1,..., « and f(x) be continuous functions in an inter- 
val ]a,b[, with a n (x)^ 0 in that interval, and let x 0 be arbitrary real values. 

Then there exists a unique solution y = y(x) of the differential equation (1) in all the interval ]a,b[ 
verifying: 


t(^o) = To ? t ? ( x o) = Ti?- • t” _ 1) ( x 0 ) = y r 
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Definition 2.2 An homogeneous linear differential equation of the nth order is a differential equation 
(1) where f(x) = 0. 1.e.: 

a „ (x)/ ) (x) + a„_ 1 ( x )/“ 1) (x) +...+ a i (x)y'(x) + a 0 (x)y(x) = 0 , (3) 

which can also be written like this: 

Ay = 0 , A = ^ j a i (x)D‘ . (4) 

i = 0 

As a consequence of the linearity of the derivative operator, the following result is obtained: 

Theorem 2.2 If {y x (x) ,y 2 (x),...,y^ (x)} are solutions of (3), then any linear combination of these 
functions of the form: 

y = C y x (x) + C 2 y 2 (x) + ...+C n y n (x) , C,,C 2 , ... ,C n eR, 
is a solution of (3). 

Theorem 2.2 suggests the following definition: 

Definition 2.3 A set of functions {y x (x),y 2 (x),...,y n (x)} , all of them continuous in the interval 
is said to be linearly dependent in ]a,b[ if there exist constants a 1 ,a 2 ,...,a n , not all naught , 

so that 

«1 y, (*) + a 2 y 2 (x)+ ... + a n y n (x)=0 , 

in all the interval Otherwise the functions { y , (x),y 2 (x),. . y n ( x ) j are said to be linearly 

independent. 

Definition 2.3 and theorem 2.2 enable the obtaining of the following result: 

Theorem 2.3 If {_>>, (x) ->y 2 (x),... ? y n (x)} are n functions linearly independent in the interval ]a,b[, 
solutions of (3), where all the functions a t (x) , i = 1,2,..., n are continuous in the interval ]a,b[ with 
a n (x)^ 0 in the mentioned interval , then the general solution of equation(3) is given by: 

y-C x y x (x) + C 2 y 2 (x) + • • • + C n (x) , 
where C x , C 2 , . . . , C n are arbitrary constants. 

Let’s observe the importance presented by the linearly independent sets of solutions of equation (3), 
{Li ( x )?T 2 (*)>• * yn ( x )} • This brings about that these sets are given a special name. 

Definition 2.4 The n solutions {y x (x),y 2 (x),...,y n (x)} to (3) are said to form a fundamental sys- 
tem of solutions of equation (3) if they are linearly independent over the definition interval of differ- 
ential equation (3). 

The following concept will enable us to characterise the fundamental sets of definitions in a simple 
way, without having to resource to definition 2.3. 
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Definition 2.5 Let jj, (x),y 2 (x),. . v„ (x)j be n functions having derivatives up to order n-1 in the 
interval \i,b\. The Wronskian of { v, (x), v 2 (x),...,_y n (x)| is defined, evaluated in x, and the value 
of the following determiner is denoted as W (y l ,y 2 ,...,y n ;x^ : 

y^x) y 2 (x) ... y n (x) 

y\(x) y\(x) ... y' n (x) , 

y n ~ l \(x) y n -\{x) - y n ~\{x) 

The way of characterising fundamental sets of solutions is given by the following result: 


w (y 1 ,y 2 ,...,y n ;x) = 


Theorem 2.4 The n solutions {>’, (xj,y 2 (x),...,y n (x)j to (3), defined in the interval ]«,/?[ are a 
linearly dependent set, if and only if 

W(y,,j 2 ,...,y„;x) = 0. 

Furthermore, if the n solutions {y l (x),y 2 (x),...,y n (x)} to (3) are continuous in the interval 
and the function a n (x) ^ 0 , then only one of the following possibilities is verified: 


1 • W(_y,,_y 2 ,...,_y B ;x) = 0 , \/x&]a,b[, 

2- W(y l ,y 2 ,...,y n ;x)*0 , Vxe]a,b[. 

We will call homogeneous linear differential equation of the n th order a differential equation (1) 
where: 

a. [x) = a i , a t eR , i = l,2,...,n ,/(x) = 0, 
and which is written in this way: 

a n y n) ( x ) + a n- 1 y n ~ l) (x) +...+ a 1 y'(x) + a 0 y(x) = 0 . (6) 

By using the derivative operator D , equation (6) results like this: 

P(D)y = 0, P(D) = a n D n) + a n _ x D n ~ l) + . . . + a x D x) + a 0 D , (7) 

For equation (6), using theorem 2.1 the existence of solutions is guaranteed. Furthermore, the follow- 
ing result will enable us to find solutions of (6) which have an exponential form. 

Proposition 2.1 The differential equation P(D)y = 0 has function y(x)=e rQX as solution, for each 
one of the roots r Q of the polynomial equation: 

P(r ) = 0. (8) 

(8) is said to be a characteristic equation of (6). 

Let us distinguish different cases according to whether the roots of (8) are simple or multiple. 
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Theorem 2.5 If (8) has n roots different form each other, real or complex, r l9 r 29 ...,r n , then the func- 
tions 

{ n x r 9 x r„ x | 

,6? 2 J, 

form a fundamental set of solutions of (6), and therefore the general solution of equation (6) is given 
by: 

y = C x e riX +C 2 e r2X + ... + C n e r * x , 
being C x , C 2 , . . . , C n arbitrary constants. 

Note 2.1 Let us observe that if (8) has a complex root a + i/3 , it also has its conjugate a -if as a 
root. Therefore, in the general solution of (6) obtained in theorem 2.5 a term of the following form 
must appear: 

Cj e ( “ +i/?)x + C 2 e ( ““ i/?)x , (9) 

with C, , C 2 arbitrary constants. If we are interested in the real solutions, we can used the identity 
e («+ib) _ Q a ^ cos ^ + i s in(/?)) , (10) 

to write (9) as follows: 

C, e ( “ +i/?)x + C 2 e ( “- i/?)x = e“ x (A cos(^ x) + B sin {ft x)) (1 1) 

where A = C X + C 2 and B = i(C x - C 2 ). 


Theorem 2.6 If T — r Q is a root of (8) with m order of multiplicity, then m solutions linearly inde- 
pendent associated to r Q exist, 

[e r ° x , xe r ° x , x 2 e r ° x m -'e r ° x }. 

Note 2.2 If (8) has a complex number a + ifl with m multiplicity as a root, then a -if is also a root 

with m multiplicity. Applying theorem 2.6 to the 2m solutions corresponding to these two roots and 
working as stated in note 2. 1, it can be proved that 

\e ax cos(p x),xe ax cos (/? x), x 2 e ax cos(/? x),...,x m ~ l e ax cos(/? x), 

e a x sin (/? x),xe ax sm(p x),x 2 e ax sin (fi x),...,x m ~ l e ax sin (fi x)j , (12) 

is a set of 2m independent solutions of differential equation (6). 

3. Solution of the homogeneous linear differential equation of « th order 

with constant coefficients with DERIVE. Implementation. Examples. CONST.MTH 

KRONECKER_VECTOR(s) := VECTOR (KRONECKERfs , 0), i, 1, DIM(s)) 

i 

VECTOR_AUXrLrAR_lCs J i) := 

If i = 1 
1 

If sfi t 0 

If sf(i - 1) t 0 
1 
0 


0 


D 
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DIM(v) - k + 1 

AUX_0(v) := VECTOR(x , k, 2, DIM(v) + 1) 

AUX_l(v) := v.AUX_0(v) 

AUX_3 (v) := SOLUTIONS CAUX_l(v), x) 

SEPARAJtAIZ_l(v) := VECTOR (IF (IMCCAUX_3Cv)) ) A 0, (AUX_3(v)) , 0), i, DIM(AUX_3 (v))) 

i i 


SEPARA_RAIZ_2Cv) := ZCKRONECKER_VECTORCSEPARA_RAIZ_lCv))) 

fd 


AUX_4(v, i) := VECTOR 


Tim 

,x^(AUX_3 (v))j,i 


(tJ 1 


AUX_l(v), j, 0, DIM(v) 


-4 


AUX_5_l(v, i) := KR0NECKER_VECT0RCAUX_4Cv, i)) 
AUX_5 (v, i) := VECTOR_AUXILIAR_2(AUX_5_lO, i)) 


AUX_REAL_l(v, i) := VECTOR 


(AUX_3(v)) -x 
k - 1 i 

c -x -e , k, 1, DIM(v) 

l i ,k 


AUX_REAI 2(v, i) := AUX_REAL_l(v, i).AUX_5(v, i) 

SEPARA_RAIZ_2(v) 

AUX_6(v) := I AUX_REAL_2fv, i) 

i=l 



r r (aux_3(v)) .x'l 

k - 1 i 

1 

AUX_CO MP L E 3 A_1 ( v , i) := VECTOR 

d -x -ReLs ^ 

, k, 1, DIM(v) 1 


AUX_CO MP L E 3 A_2 ( v , i) := VECTOR 


i >k 


■ SIGNClMCCAUX_3(v)) )) ■ IM 1 


f (AUX_3 (v)) .xT 

L f J. 


k, 1, DIM(v) 


AUX_CO MP L E 3 A_3 ( v , i) := AUX_COMPLE3A_l(v, i) + AUX_COMPLE3A_2(v, i) 

AUX_COMPLE3A_4(v, i) := AUX_COMPLE3yL3 (v, i)-AUX_5(v, i) 

AUX_7(v) := I(VECT0R(AUX_C0MPLE3A_4(v, i), i, SEPARA_RAIZ_2(v) +1, DIM(AUX_3 (v)) , 2)) 
DTFERENCIAL_HOMOGENEA_COEFICIENTES_CONSTANTE5(v, x) := EXPAND(AUX_6(v) + AUX_7(v)) 

fd 'ij 

V_C_AUX_2fv, c, i, ]) := I — I c 

V_UVUX_3Cv, c) VECTOR^— J"* c, j, 0, DIM(v) - 2 


-1 

V_C_AUX_4(v, c) := V_C_AUX_3(v, c) 


V_OVUX_5Cv, c, f) := VECTOR 


f 

KRONECKERCj, DIM(v) - 2), j, 0, DIM(v) - 2 


v 


1 


V_C_AUX_6(v, c, f) := V_C_AUX_4(v, c)-V_C _AUX_5(v, c, f) 
V_UVUX_7Cv, c, f) J V_Cj“iUX_6(v, c, f) dx 
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V_CJWX_S(v, c, f) := VECT0RC1, j, 0, DIM(v) - 2) 
V_C_AUX_9(v, c, f) := V_C_AUX_7(v, c, f).c 


SO L U CIO N_P ARTI CU L AR_C0 MP L ETA ( v , c, f) := V_C_AUX_9(v, c, f) 


Given the homogeneous differential equation of the //th order with constant coefficients 

a n y n) { x ) + a n-i y n ~ l) { x ) + ••• + ay(x) + a 0 y(x) = 0, 

we will introduce in the vector v = [a n ,a n _ l9 ...,a 0 ] the coefficients of this equation, written in serial 
order, tapering downward. The CONST.MTH file functions operate in the following way: 

1 . Function A UX_0(v) creates vector \jc n , x n ~ l , . . . , jc, 1 ] with the same dimension as v. 

2. Function A UX_\ (v) starting from v and^C/X_0(y) provides us with the characteristic polynomial of 
equation (6). 

3. Function AUX_2(v ) calculates the roots of the characteristic polynomial (let us remember that De- 
rive does not directly provide us with the multiplicity of a roots of a polynomial). Thus AUX_ 2(v) 
gives a vector of the form [x = r x , x = r 2 , . . . , x = r s ] , where [r x , r 2 , . . . , r s } are the roots of the characteris- 

s 

tic polynomial, real or complex, each one with a multiplicity {/ 1? / 2 ,...,/ 5 } and being verified ^J, i = n. 

i = 1 

4. Function AUX_ 3(v) provides us with the list of roots calculated by the function A UX_2(v), giving a 
vector of the form [r 1? r 2 ,...,rj . 

5. Functions SEPARA RAIZ l (v) and SEPARA_RAIZ_2( v) separate the real and complex roots of the 
characteristic polynomial given by AUX_3(v). 


6. Function AUX_4(v,i) , for every root of the characteristic polynomial r , give us a vector with the 
values of the characteristic polynomial and its n - 1 first derivatives of that root. Thus it produces a 
vector of the form [b 0 , b x , . . . , b n _ x ] , where 

d J P(x) | 


b.=- 
1 dx J 


, 7=0,1,. ..,72-1. 


It is evident that when a root has multiplicity 5, it will give rise to a vector of the form 


0,. . .,0,Z^ +2 ,. . .,b n _j 


7. Starting from a vector 


0, . . . , 0 , b s+2 , . . . , b n _j 


, obtained by means of the function AUX_ 4(v,/), the 


function A UX_5_l(v,i) provides us with a vector of the form: 


s+\ 

0,0,.. .,0 


for each root I) of the characteristic polynomial. For that purpose it uses functions 
KRONECKER VECTOR (s), VECTORA UXILIAR_\(s,i) and VECTOR_AUXILIAR_2(s,i). 
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8. Function AUX_5(vJ ) filters function AUX_ 5_l(v,i) so as to avoid anomalies such as the one in the 
following example: 


Example 3.1 Let us consider the polynomial P(x) = x 4 + x 2 . This polynomial verifies 


- x 4 + x 2 


4x + 2x 

2 


P(x) 

r(x) 

P"(x) = I2x 2 +2 
P”'(x) = 24x 
P 4) (x) = 24 


P(0) = 0 
P'(0) = 0 
P"( 0) = 2 
P"(0) = 0 
P 4) (0) = 24 


■ x = 0 is a double root of P(x) , 


Function AUX_4(v,i) for v = [1,0, 1,0,0] provides as a result vector [0,0,2,0,24] and function 
AUX_5_l(v,i) provides vector [l,l,0,l,0]. Thus strange solutions would be introduced to differential 
equation (6). Function AUX 5 (v,i) permits us, for this concrete example , to obtain vector [1,1, 0,0,0] 
as an answer, which interprets the root x = 0 as a double root of this polynomial. 


9. Functions AUX_REAL_l(y,i) and AUX_REAL_2(v,i) work with the real roots, operating with each 
one of them so that they appear with its corresponding multiplicity. 

10. Function AUX_6(v,i) gives us the partial solution of the differential equation obtained of the real 
roots. 

11. Functions AUX_COMPLEJA_\(y,i), AUX_COMPLEJA_2(v,i), AUX_COMPLEJA_3(v,i) and 
AUX_COMPLEJA_4(v,i) operate on the complex roots r z of the characteristic polynomial. Basically 
they operate their real and imaginary parts separately, and prepare them so that when they are added in 
pairs, as indicated in notes 2.1 and 2.2, the imaginary numbers disappear and the root appears with its 
corresponding multiplicity. 

12. Function AUX_l{y) adds all the partial results obtained when working with the real and complex 
roots separately, providing us with the solution of the differential equation. 


Examples: 

y'”(x) + 3y”(x ) -4/(x)-12j(x) = 0 

In this case vector v = [l,3,— 4, — 12] , and we have: 

DIFERENCIAI H0M0GENEA_C0EFICIENTES_C0NSTANTES([1, B, -4, -12], x) 

2 ■ x - 2 ■ x - 3 ■ x 

e -c + e •€ + e -c 

1,1 2,1 3,1 


Let us check that the solution obtained is actually a solution of (13). 

2 ■ x - 2 ■ x - 3 ■ x 

y(x) := e « c + e ■ c + e - c 

1,1 2,1 3,1 


"d " 

3 

"d ” 

dxj 

y(x) + 3- 



d 

y(x) - 4- — y(x) - 12 ■ y(x) = 0 
dx 


(13) 
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The next two examples will be solved by applying CONST. MTH and by MATHEMATICA: 

y (4) - 7/-18y = 0 
y (4) +y" = 0 

DrFERENCrAL_HOMOGENEA_COEFICIENTES_CONSTANTESC[l, 0, -7, 0, -18], x) 

3 ■ x - 3 ■ x 

e -c + e *c + C0SC/2-x)-d + SING/2 ,x),r 

1,1 2,1 3,1 3,1 

DI F E R E N Cl AL_H0 MO GE N E A^CO E F I Cl E NTE S_C0 N STANTE S ( [ 1 , 0, 1, 0, 0], x) 

C0S(x)-d + SIN(x),r + x-c + c 

2,1 2,1 1,2 1,1 

in [i ] := DSolva [y T T 1 7 [x] - 7 y 1 T [x] - 18 *y [x] >■ 0 , y , x] 

out [i ]= {(y — ► Function^ [x], •■e~ 3x C3 + C 4 + c\ cos (V2 x) + C 2 sin(V 2 x)]| | 

in[ 2 ]:= DSolva [y T T 1 1 [x] + y T T [x] == 0 , y , x] 

uijt [2]= | !'_y -s Fiuictioii[!x], c? + x C4 - c\ cos(x) - C2 sin(x)]}} 


Finally two more examples of order 5: 


y 5) + -y 4) + — y m + ^^y” + — y + — y = 0 


237 1358 „ 108 , 72 


25 


15 25 125' 


_y (5) - 3y m + 4y' = 0 

DERIVE calculation + check of the solution: 

DIFERENCIAI HOMOGENEA_COEFICIENTES_CONSTANTES 



237 

1358 

108 

72 ' 

xl 

{[ ■ s ’ 

25 

125 

25 

125 . 

' J 


2 - 2 ■ x/5 - 2. x/5 - 2 ■ x/5 

x -e -c + x-e -c + e -c + C0S(3-x)-d + SIN(3*x)*r 

1,3 1,2 1,1 2,1 2,1 

2 - 2- x/5 - 2- x/5 - 2- x/5 

y(x) := x -e ■ c + x-e -c + e -c + CQS(3-x)'d + SIN(3*x)»r 

1,3 1,2 1,1 2,1 2,1 

6 ■ y ' ' "(x) 237-y rr, (x) 1358* y"(x) lOS-y'(x) 72-y(x) 

y' ,p ''(x) + + + + + = 0 


25 


125 


25 


125 


DI F E R E N Cl AL_H0 M0 GE N E A_CG E F I Cl E NTE S_C0 N ST ANTE S ( [ 1 , 0, -3, 0, 4, 0], x) 

■ d 


y(x) e 


f X - ] yZ-x/2 f xl - y7-x/2 f xl - V7-X/2 f x ■) 

■ COS \+e -r -SIN +e *d ■ COS \+e *r ■ SIN +c 

,1 L 2 J 2,1 l 2 J 4,1 L 2 J 4,1 l 2 J 1,1 

f x ^ 7? ■ x/2 fx) - J7* x/2 C x ^ x/2 f x 1 

■ d -COS — \+e -r -SIN — +e -d -COS — +e -r -SIN — +c 

2,1 L 2 J 2,1 l 2 J 4,1 l 2 J 4,1 l 2 J 1,1 


47 -x/2 


y (x) - 3*y' p, (x) + 4-y'(x) = 0 


[1] T.M.Flett, Differential Analysis. New York. Cambridge University Press. 1980 
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THE RIEMANN INTEGRAL WITH THE TI-92 

Francisco Jose Santonja Gomez 

C.U. ESTEMA-ANTONIO DE NEBRIJA UNIVERSITY. Valencia. Spain 


In this paper some Tl-programs and functions concerning integration are given. 


Description of the Tl-programs Firstly the Lower Riemann Sum is 

shown by cajainfQ. 

: eajainfl 3 Let fix) = x 2 on [0,3] with 5 subdivi- 

: Prgn sions. 

: setHode(" Exact/Approx", "EXACT") 

: D l a l og imag 

: Ti t T e ''INPUT DATA" , Irrtsryal tourah B 

:Reque$t "Interval . 1 ower bound", a i"!"?* 1 UPf ~ r i baurtd: k 

: Request "Interval .upper bound" t b . " 

Request "Subinterval s" ,n i»l 

:Request "function" ,f 

: EndDl og M?==^ tm --irwg 

: exprta )-*e ;expr(b l^brexpr t n 3->n :expr ( f )->f n T~^r 

;C1 r i — KOQH |Tr-ace iRgfr aphrlithlo^iw 

:For 1,0, n ( l P ; ; , 

: Pt On e+( b-a 3*1 /n , 0 [ ; ; : /_ 

; : EndFor \ ■ 1/T 

: 1 tmitt f , x, m)-»val or( f ,tn) , : 

: For i,0,n,l ■ • ■ 

:Pt0n a+( b-a 3*1 In , val or ( f ,a+(b-a)*i / n } ■ • — L 

: E n d F o r mHI — mt ehso " fu*k 

: F Q T I 0 „ 11 A 1 /h. A U{/ n ] 

; EndFor 
i For 1 , 1 , n , 1 

: Li ne a+(b-a)*(i-l)/n f valortf,a+(b-a)*{i-l>/n},a + (b-a)*i/n, 
valorCf,a+(b-a)*(i-l)/n) 

: EndFor 
: EndPrgm 

Next, the program cajasup() has been defined to present the Upper Riemann Sum. 

You only have to replace the marked program block of cajainfQ by the block given below. 

^ajasupO PEEKES 

;Pr@m N 1 [ ' — s — — U ~ J — L 


kt— ifU Cal cPthgr jrrny .10 ICl_«»r a-z.. 

gjSag Z 

t Iirtsrysl lewsr bounds j£) 

Interval sjpp t r botaida I 

i5 I 

i Fwiktion* m. i * 

»[ ') tB5C=CHHCELb . 


cajaxnl 

Ihhim 




: For 1 ,0 , n-1 , 1 

: PtOn a+( b-a )*i/n, valor (f ,a+< b-a) *(1-1 )/n) 

; EndFor 

: 1 1mltff ,x .iu>H>val or( f ,m) 

.-For I.O.n-l.l 

: Li ne a+{ b-a 3*1 /n , 0 , a+( b-a 3*1 /n , 
Yalor(f,a-Mb-a>#(i-i-l>/n> 

: EndFor 

: Line b,0 , b , val orCf ,b) 

: For i , Q ,n- 1 

: Line a-M b-a )*i /n , val or ( f ,a+( b a ) + £ 1+1 ) fx\ ) 
a+(b-a)*(i+l)/n, yal or ( f, a-Hb -a )*E 1-i-mn 3 
; EndFor 



Combining ca j a 1 n f ( 3 with cajasupt ) gives 5 up 1 Tift ), which presents both sums (showing above 10 
subdivisions). 
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Next program rieman() has been defined to calculate both Riemann Sums. 

: MematU ) 

: Prgm 

: Dialog 

;Title "INPUT DATA" 

Request "Interval . lower bound", a 
: Request "Interval . upper bound" , b 
{Request “Subinterval s" ,n 
: Request "functi on" , f 
: EndDlog 

:expr(a )-»a :expr C b}*biexpr f n )->n :expr( f )-*f 
:liniittf,x T m)-*valor(f,n> 
tinax(valor(f f a> f valor(f,b) )-*hina( f , a , b) 
mini vol or(f ,a ) .valor C f.b) ( f , a , b ) 

{sum($eq((b-e)/n*(Jim1 <f , a+< b-a)*( 1 -1 )/n ,a+(b-a )*i/n ) ) , t , 1 ,n) )->f nf surma 
: setModat "Exact/ Approx" , "APPROXIMATE" ) 

!sum(seq((b-a)/n*(tima (f , a+( b-a)*l i -1 )/n,a+(b-a )*l/n ) ) , t , I , n ) Wsupnuma 
: Dialog 

{Title "RIEMAHN SUMS" 

:Text "Lower sura M &string(l nfsuna } 

;Text "Upper ng< supsuraa ) 

: EndOl og 
: EndPrgm 

Finally program tableta() shows the integration as the limit of the sums (lower and upper sum) con- 
verging to the same value. 



* = 10 


:tableta ( ) 


Prgm 

ClrTable 
Cl r 10 

$ etModeC " Exa ct/Approx" , "APPROX I HATE" ) 
Dialog 

Title "INPUT DATA* 

Request "interval . lower bound", a 
Request "interval , upper bound” ,b 
Request “function" .f 
EndDl og 

expr(a )*a ;expr( b)->b:exprt f )->f 



: Table $UTH{$eq((b-a)/n*(btni(f,a+(b-a)*(i-l>/n,a+(b-a)*i/n)),i,l ( n)> J n 
{Table surnt $eq( ( b-a 1 /n*( hmatf , a+( b-a )*( i -1 >/n , a+( b-a)*i /n) ) , i ,l,n)) ,n 
: End P rgm 


References: 


Llorens Fuster, Una Lection de Matematicas con Ordenador, EPSILON, n.31-32, p. 81-88. 1995. 
Llorens Fuster, Aplicaciones de DERIVE al Analisis Matematico 1, Polytechnic University ofValencia 
Press. 


When l received F Sanionjm Hies i immediately tried them because some times ego - at the occasion of the 
legendary Kmms Conference 1 in i§92 - l produced a DERIVE package on the same topic, {Happy memories on 
Terence Etchetls' pm email - correspondence and diskette exchanging by snait mail}. I found the TJ-programs 
useful but had the idea that they could he more userfriendiy, M&ny settings must be set from outside of the pro- 
gram: graph the function, set me [window] values., mi the [TciSet] for the last program, ta bl eta E ) does not work 
without running r 1 eman ( ) before, And look what is happening when using a non increasing function: 



These packages could sen/e excellent as & starting point for stu- 
dents who rant to improve their idem to represent and to calculate 
Riemann sums. Step by step they could try to Include additional 
features to the core of the programs and at last try to combine 

ihem„ 







For me this contribution was challenge enough to remember my DERIVE deeds from 1992 and I pro- 
duced a TI-92 package: :integ(). My friend Wolfgang Propper saw my efforts and very similar to me he 
again found many ways to make this package more user friendly and he added some important fea- 
tures. Encouraged by Bernhard Kutzler we produced a book "The Riemann Integral - From Counting 
Raindrops to the Fundamental Theorem". I’d like to add some screen shots to show where F. San- 
tonja’s inspiring ideas led Wolfgang and Josef: 





Shew the graphic representation of the methods: 
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30 
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m 
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3,3314 

5.25SG 

M 


5.4139 
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S.3&4* 

S*2S04 

UMlI 



Compare the convergence by inspection: 



Work symbolically: 



Discuss the integral functions: 


(... and much more.) 
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Helge von Koch' s Snowflake 


Matija Lokar, Faculty of mathematics and physics, University of Ljubljana, Slovenia 

Matiia.Lokar@fmf.uni-lj.si 


As at the time being it is snowing, let’s try to draw some snowflakes with DERIVE, too. As basis we 
take the so called Koch curve. It belongs to the class of recursive defined objects. We get them with 
applying a certain procedure to a base object. Base object, which can be line, box, circle, triangle or 
something else, is transformed into several similar objects which serve as the starting point of thea 
same transformation. We know Peano curve, Dragon curve, Hilbert curve, Sierpinski triangle and 
several others. 



We start with line (Koch curve of order 0). The basis procedure is as follows. We divide the line into 
three parts. The middle part is substituted with two lines that form en equilateral triangle with deleted 
line. We get the Koch curve of the 1 st order. If we apply the procedure on all four segments then we 
will get the Koch curve of 2 nd order. 

An interesting paper about Koch curve can be found on Internet on page 

http : //www . euro logo . org/papers/logomov . html [1} . 


There LOGO is used. But we would like to draw with DERIVE. Let us first recall some basis facts 
about drawing with DERIVE. 

A point in DERIVE is represented as pair of coordinates in square brackets, divided by a comma. So 
the origin is represented as [0,0] and a point with co ordinate x = 1 and coordinate y = 2 as [1,2]. If we 
plot such an expression ([1,2] for example), we get point (1,2). If we would like to draw several 
points, we make a vector consisting of two-dimensional vectors, namely an n x 2 matrix, where n is 
number of points. When we draw several points, DERIVE can draw a line which connects points. 
When we are in 2D-Plot Windows, we choose Options State and set Mode to Connected or 
Discrete. We can also set the size of points. When set to Small, point can not be distinguished 
from the joining line in Connected mode. 

Now to our Koch curve. First we calculate three new points that our procedure makes. Each one will 
be calculated with a function. The parameters for all three functions are the starting and ending point 
of the line. We exploit the fact that we can calculate with points like with vectors. 

TT (a,b) : =a+ (b-a) /3 
TDT (a,b) : =TT (b, a) 

TO (a,b) := (a+b) /2 + SQRT (3) /6* [a SUB 2-b SUB 2,b SUB 1-a SUB 1] 


Let us check first our attempt. Set u as a starting point ([-1,-1]) and 
v as an ending point ([1,1]). We have to draw a polygon line 
[u, TT (u, v) , TO (u, v) ,TDT (u, V) , v] . 

It seems OK. Now we use the recursion. If the order of curve is 0, 
then we have a line. Otherwise we have four curves of one less 
degree on: u-TT (u, v) ; TT (u, v) -TO (u, v) ; 

TO (u, v) -TDT (u, v) and finally TDT (u, v) -v. 



D-N-L#33 


Matija Lokar: Helge von Koch’s Snowflake 


P 47 


We edit and plot KOCH ([-2,0], [2,0], 3]. Simplifying we 
get a huge expression. Before plotting we turn off automatic color 
cycling and hide the axes (all in Options Display). 



But if we observe the received expression, we are not quite satisfied. First there is no need to use exact 
arithmetic, so we set Precision := Approximate PrecisionDigits := 2. This also 
significantly speeds up the computation. Also on each step we end up with four polygon lines. But we 
should get only one. So we change function KOCH a little: 

K0CH(a, b, n) := If n = 0 
[a,b] 

APPEND(KOCH(a,TT(a, b) , n-1) , K0CH(TT(a, b) ,T0(a, b) , n-1) , 
K0CH(T0(a,b) ,TDT(a, b) , n-1) , KOCH (TDT(a, b) , b , n-1)) 

Another improvement is due to observation that we get all points twice - first as an ending point and 
secondly as a starting point. If we just drop the second point in case of order 0 then we circumvent this 
obstacle: 

K0CH(a, b, n) := If n = 0 
[a] 

APPEND(KOCH(a,TT(a, b) , n-1) , K0CH(TT(a, b) ,T0(a, b) , n-1) , 
K0CH(T0(a, b) ,TDT(a, b) , n-1) , KOCH (TDT(a, b) , b , n-1)) 

Now our function does not behave “properly” for n = 0 and we lost the first point of the sequence. So 
we define: 

KOCH_CURVE(a, b, n) := IF(n=0 , [a, b] , APPEND ([a] , K0CH(a, b, n))) 

If we draw three Koch curves based on the sides of an equilateral triangle the we get the so called 
Koch Snowflake . First let us define POLY (n) which produces the vertices of a regular polygon of 
order n. 

POLY(n) := VECT0R([3-C0S(i) ,3-SIN(i)] ,i , 2 tt, 0 ,-2/r/n) 

Then we define KS (a, m) which creates a collection of KOCH curves of order m between consecu- 
tive points of a vector a. 

KS(a, m) := VECTOR (KOCH (a^i ,aj,(i+l) , m) ,i ,1, DIMENSION (a) -1) 

Simplifying KS (POLY (3) ,5) gives the Snowflake. 

We see it together with an accompanying triangle. 

It should be also noted that the order n in which we 
compute the polygon is important. If we use points 
computed with 

POLYB(n) := 

VECT0R([3C0S(i ),3SIN(i)],i , 0 , 2 *7r, 2tt/ n) 
then we get Koch curves which are inside the triangle. 
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KS(P0LYB(3) , 4) 


KS(P0LY(6) , 3) 


KS(P0LYB(6) , 3) 




The snowflake is always a challenge to be faced. I think that it may be nice to compare the 
various ways producing this famous fractal structure. You can do this in this DNL studying 
the next contribution. Compare also with the Roanes-Lechner-Wiesenbauer implementation 
of LOGO in DERIVE (in DNL#25, March 1997) and with Robert Setif s Treasure Box (DNL#3 
and DNL#6 from 1991/92). Josef 


Another comment from 2014 concerning “challenge”: I tried to reproduce the recursive pro- 
cedure given above with TI-NspireCAS. The first steps - were not so difficult - even missing 
the programming tools for graphics included in the V200 programming language. 

But I could not find a way to transfer Matija’s recursive creation of the Snowflake. It would be 
great if any Nspire expert among you would accept the challenge and follow Matija’s ideas 
from 1999 with the nowadays TI-NspireCAS, Josef 

This is what I did so far: 


Define koch(a, b,n)= 

Prgm 
Local k 
If n=0 Then 

kl : =col Au gmentfa, b ) 

Else 

kl : =ko chfcol Au gmentffr 7, tt{a, 7?) ), n - 1 ) 
kl : =ko chfcol Au gmentffr 7, to {a, £?))//- 1) 
kl : =ko chfcol Au gmentf Jc7, ttftia, b) ), n - 1 ) 
kl : =ko chfcol Au gmentffr 7, b) f n- 1 ) 

End If 
kl 

End Prgm 



[1} Now this website can be accessed at https://eurologo.web.elte.hu/lectures/logomov.htm 
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Computergrafik mit DERIVE - Computer Graphics with DERIVE 

Dr. Maria Koth, University of Vienna, Austria 


2.2 All diagonals of a regular n - gon 

A nice graph is created connecting each vertex of an 7z-gon with each of the remaining vertices: 

The radius of the circum circle is 1. The line connecting the / lh vertex [COS (i) , SIN (i) ] with the 
j th one [ COS ( j ) , S IN ( j ) ] is described in DERIVE by a list 

[ [COS (i) ,SIN(i) ] , [COS (j) ,SIN(j) ] ] 



To avoid plotting lines twice it is sufficient to 
connect the / th vertex with vertices # (7+1), (7+2), 
....n. Function LINIEN ( i , n) creates all the 
lines starting from point i. 

DIAGON (n) produces all diagonals by running i 
through all n vertices of the n- gon. 


You can see the result of DIAGON (20) . 


LINIEN (i, n) := VECTOR ( [ [COS (2*7T*i/n) , SIN (2*7T*i/ n) ] , [COS (j) ,SIN(j) ] ] , j , 

2*tt* (i+1) /n, 2*71, 2*7i/n) 

DIAGON (n) : ^VECTOR (LINIEN (i, n) , i,n-l) 

2.3 A family of rotated n-gons 

These graphs show families of polygons which are produced by rotating a polygon by a constant angle 
8. [ eg. 8=5° and simultaneously shrinking it in such a way that the new polygon is exactly inscribed 
the preceding one. At first one has to find the shrink factor /: 

Given is sector OPiP 2 of a regular 72 -gon with radius r. Its angle a in the 
center is a — ZP 1 OP 2 = lain. Angle 6 is given by fi = (180° - a)/2 = 

= tt/2 - nln. The rotated triangle ORiR 2 is part of the next polygon with 
a circumscribed circle of radius f-r = ORi. We apply the sine rule in 
order to find ORd 

rf r sin fi sin( k! 2- n / n) 

sin fi sin (1 80 °-j3 - 8) sin(l 80° -ft -8) sin(;r 12 + 7t / n - 8) 



Using sin(90° - a) = cos a gives /: 
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F(n s 5):= 



COS 



Setting r = 1 gives the radius of the 
k th polygon’s circumscribed circle 
with/. This polygon is turned by the 
angle k-S. 


ROTECK returns its vertices. 

ROTS CHAR (n, 5 , k) produces a family of z rotated n-gons. 

ROTECK (n,8, k) := VECTOR ( [F(n, §) A k-COS (i + k-8) , F(n,8) A k • SIN ( i + k • 5) ] , 
i, 0, 2 *71 , 2 *7i /n) 

ROTSCHAR (n, 5, z) == VECTOR (ROTECK (n, 5, k) , k, 0, z) 

ROTSCHAR ( 4 , 5% 20) 

ROTSCHAR (10 , 30° , 20) 

ROTSCHAR ( 3 , 2 • ° , 25) 

3. The Koch Curve 

The Koch curve is a well known example of a fractal. It was defined in 1904 by the Danish mathema- 
tician Helge von Koch. Koch had the intention to show a ’’Monster curve”, which is continuous on 
each place but nowhere differentiable. 

3.1 The construction 

We assume that you all know how to obtain this famous curve: 

Using DERIVE commands you can produce a 
Koch curve or order n in an easy way. You start 
with the unit segment as Koch curve of order 0 
and build step by step curves of next order. The 
important idea is to relate a segment (x,y) to a 
sequence of segments generator (x, y) . 

The only difficulty is to express the vertex of 
the equilateral triangle by the points x andy. We 
define the auxiliary function n v ( v ) , which is 
the normal vector of v with equal length and 
turned to the left. Using nv ( v) we can define the generator - function. 



D 
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nv(v) := v- 


0 1 

-1 0 


generatorfx , y) := 


2- x + y x + y 


x 3 


next(v) := APPEND(VEG"OR(generator(v a v 

i i 


(x+y)/2+NV(y-x)sqrt(3)/6 



x 


y 


J 3 x + Z ■ y 

+ -nv(y - x), 

6 3 


) a i i DIMCv) - 1), [Tv 11) 

+ 1 LL DIM(v) J J 


koch(l) := 

If 1 = 0 

[0 S 0; 1, 0] 

APPEND(nextCkoch(l - 1))) 


The Koch curve of order n consists of a list of points v = [vi, v 2 , ... v D im(v)]* You receive the curve of 
next order n + 1 in applying generator to each of the segments (vi. Vi+i). next (v) creates these 
vectors generator (v if v i+ i) and appends them one after the other. Finally we have to append the 
last point v D im(y). Starting with segment [[0,0], [1,0]] as Koch curve of order 0 we use a recur- 
sive construction process to create Koch curves of order k with koch ( k) . 

I will demonstrate how generator works with the next curve (Koch2-Curve), Josef. 


Even if the application of matrices to describe mappings in the 
plane has not been discussed in classes you could provide this 
function as a "Black Box". 


koch(l) 
koch(Z) 
koch (3) 


koch (4) 


3.2. Koch similar curves 

You can generalize the Koch curve in various ways. If you replace generator (x, y) by any other 
appropriate function, koch (k) will generate "Koch similar” curves. Four examples of these curves 
will be presented in the following, nv and next will remain the same. 


3.2.1 The Koch-2-Curve 

If you start with the unit segment and choose points [0.4, 0.2] and [0.6, -0.2] as vertices of the generator 

three segments of equal length — . 

v 5 

This generator is described by DERIVE in the following function: 

generator (x,y) :=[x,0.6*x+0.4*y+0. 2*nv (y-x) ,0.4*x+0.6*y-0. 2*nv (y-x) ] 
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Demonstration of the “generator”: 


Let’s assume that the initial segment is [(1,1), (5,3)] - which is then koch (0) . According to 
the algorithm given above koch ( l ) =append (next ( koch ( 0 ) ) ) . 


next (koch (0) ) =append (generator (koch (0) i, koch (0) 2 ) , koch (0) 2 ) 

generator ( (l,l),(5,3)) = ((l,l),0.6*(l,l)+0.4* (5, 3) +0 . 2*nv (4, 2) , 

0.4* (1, 1) +0 . 6* (5, 3) -0 .2*nv(4, 2) 


0 . 2*(4,2)* 


v-1 Oy 


=(-0.4,0. 8) 


generator ( (1, 1) , (5, 3) ) = ( (1, 1) , ( 0 . 6+2-0 . 4 , 0 . 6 + 1 . 2 + 0 . 8 ) , 

(0. 4 + 3 + 0. 4, 0.4 + 1. 8-0. 8) ) = 

= ( (1,1) , (2.2,2. 6) , (3.8, 1.4) ) 

next (koch (0) ) = ((!,!), (2. 2, 2. 6), (3. 8, 1,4), (5,3)) 


APPEND next 


We check the p&p-calculation: 

r i i 

1 1 T | 2.2 2.6 

5 3 JJ 3.8 1.4 

_ 5 3 

The graph of the first steps of this curve is 
given on the next page. 


3 


2 


1 



1 2 3 4 5 


3.2.2 The Cross Stitch Curve 

Another possibility is to replace the equilateral triangle by a square. Then the generator is built by 5 
segments of length 1/3. 


generator (x,y) := [x, (2*x+y)/3, (2*x+y+nv (y-x) ) /3, (x+2*y+nv (y-x) )/3, (x+2*y)/3] 


3.2.3 The Ice Curve 

A Koch curve which makes us think on ice needles and ice flowers can be generated very easily. 
Let us erect a segment of length 1/3 in the mid point of the unit segment. 

generator (x, y) :=[x, (x+y)/2, (x+y) /2+nv (y-x) /3 , (x+y)/2] 


The pictures on the next page show these three varieties of Koch curves: {DERIVE 6): 


Koch-2-Curve 
Cross Stitch Curve 
Ice Curve 
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[£*1 Algebra 1 


# 11 : 


generator(x, y) := 


x 


x + 2«y + nv(y - x) 


3 


E02D-plotl:l r^l~Bira1 


koch(l) 


0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 


2»x + y 2»x + y + nv(y - x) 

9 

3 3 

x + 2«y 

3 

0.6 

0.5 

04 koch(2) 

0.3 

0.2 

o.i 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 
- 0.1 

- 0.2 

-0.3 


ET 




Ifel Algebra 1 

#12: generator(x, y) := 


x + y 


x + y 


x, 


x + y 


E02D-plotl:l 


ram® 


koch(l) 


l 

0.8 

0.6 

0.4 

0.2 


0.2 0.4 0.6 0.8 1 1.2 1.4 

- 0.2 


[J] 2D-plotl:2 


koch(3) 


0.6 

0.5 

0.4 

0.3 

0.2 

0.1 


0.1 0.2 0.3 0.' 
- 0.1 


- 0.2 

-0.3 


nv(y - x) 

> 

3 


P 

1 ° || El Ifsl . 
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3.2.4 The Cesaro Curve 

The equilateral triangle in the original Koch curve is replaced by an isosceles triangle. The length of 
its sides is calculated in order to obtain a generator consisting of four segments of equal length. Vary- 
ing the base angle a results in varieties of this kind of Koch curve. 

1 1 

r = r cos a — > r = 

2 2 (1 + cos a) 


The DERIVE - Code: 

[a: =85°, r:=l/ (2* (1+COS (a) ) ) ] 

generator (x,y) : = [x, x+r* (y-x) , (x+y)/2+r*SIN(a) *nv (y-x) , y-r* (y-x) ] 
KOCH (6) 



3.3. Koch curves in a square 

Four Snowflake curves can be assembled to form a ’’Koch Square”: 
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To generate a Koch- or Koch similar curve on any chosen segment v needs only a change in 
koch ( k) from above: 

kochl ( k, v) : =IF ( k=0 , v, APPEND (next ( kochl ( k-1 , v) ) ) ) 

Using the vector v 0 to describe a square, square ( k , vO ) creates the square of order k. 

v0: = [ [0,0] , [1,0] , [1,1] , [0,1] , [0,0] ] 

square (k,v_) : =VECTOR ( kochl (k, [v_SUB i,v_SUB (i+1) ] ) ,i,DIM(v_) -1) 

Depending on the chosen generator you can receive the following pictures. Using v\ returns curves 
which are directed out of the square. Interesting pictures result in combining both versions. 

vl: = [ [0,0] , [0,1] , [1,1] , [1,0] , [0,0] ] 

generator (x,y) := [x, (2*x+y)/3, (x+y) /2+SQRT (3) /6*nv (y-x) , (x+2*y) /3] 
square ( 4, vO) vl: = [[0,0], [0,1], [1,1], [1,0], [0,0]] 


square ( 4 , vl ) 



(Students could try to find their own generators. See here some suggestions. Josef) 
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What about the Koch curves and TI-NspireCAS? 


lines 

Define lines!/, tf)= 

Prgm 

MD) 

2' 31' (j+ll 2' JT 

For ], 2 ■ jt, 

n n 


Is: 


=augment|/s J | cosj " ^ * j^inj " ^ * j,cos j)),sin(/) Jj 


End For 
End Prgm 



475 


diagonals 


Define diagonaJs(^)= 

Prgm 
Local d 
<*:={[]} 

For i, 1,^-1, 1 
lines (/rc) 

ds: =augment(tf5,/5) 

End For 
xp : =seq (cfr [i ], / 1, dim (ds), 2 
yp : =seq (tfs [ i ], / 2, dim (cfr), 2 
End Prgm 


El 


Diagonals of an n-gon 

Change the number of vertices: 


diagonals (20) ► Done 

H 
: : 


7/8 

1 


There were no problems presenting the n-gons, but things changed treating the Koch curves. 


nv(v):=\ 


■ Done 


0 1 
-1 0 _ 

ca (a ,£):=col Augment^,/?) ^ Done 

( \ |f( 2'x+y\x+y E / \\ 

generator : =ca ca ca ^ L + nv^-jm, 

i \ \ 3 / 2 6 / 


x+2 -y 


* Done 


genera tor ([ 0 0 ] J [ 1 0 ]) 


0 0 

1 

- 0 

3 

1 1L 

2 6 
2 

- 0 

3 


I tried to reproduce Maria’s recursive procedures 
with Tl-Nspire. It is a shame, I failed. It would be 
great if somebody from the DUG community could 
transfer Matija’a and/or Maria’s recursive proced- 
ures on the Nspire. 


next 

Define next(v)= 

Prgm 

m : =getierator (v[l ] , v[ 2 ] ) 

For i,2,dim(v)[l]-l 

m\=ca{m, generator (v[ i],v[ rf l])) 

EndFor 

m :=ca(m, v[ dim (v) [ 1 ] ]) 

EndPrgm 


Then I changed my strategy and used the generators iteratively - and this lucked. 
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Computer Graphics performed with TI-NspireCAS 


D-N-L#33 


I must use colAugment for collecting the points as a whole and to add y as last element of the 
generator. 


nv 


(v):=v 


■ Done ca(a^):=colAugment(fl,£) ► Done 


Adapt the generator according to your ideas: 
This gives the Snowflake: 




x+y \ _\+y+nv(y-x) \ -X+y \ 

3 \ 2 ]' 3 j 


* Done 


koch([o 0 ]j [ 1 0 ]j 4) ► Done dim(xp) ^ 626 

n > 4 is not possible because of restricted resources!! 

Lists xp and yp form the scatter diagram. 


Just for fun: what will this look like? Exchange generator and koch given above! 

generator^yJ^catcafcafcafx^x+yJ/BMx+y+nv^-xJJ^Mx+yJ/S^y) 

koch([0,0],[1,0],4) 


koch 


Define kocli(zAvyr)= 
Prgm 

Local nn,g 
fr\ =cfl(u,v) 
fro\=ft' 

If w =0 Then 
Goto end 
End If 


For i, 1 ,n 

ff/r=dim(tf"o)|_lJ 

fr=fro[l\ 


For/ l f nn-l 

g : =gene rator {fro [7 ]fro [7 + 1 ]) 
fr':=ca(fr,g) 

End For 


fro\=fr 
End For 
Lbl end 

xp : =mat ► listj/ft? T [ 1 ]) :>p : =rnat ► listffro T [ 2 ] ) 


End Prgm 




This is the generator of the fractal on 
the right hand side. 


D-N-L#33 


ACDC: ANIMATION IN DERIVE 
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Matija Lokar, University of Ljubljana, Slovenia 

ANIMATION IN DERIVE 

Motivation for this contribution is twofold: During systematic study of all menus in DERIVE 
(DOS version), I wondered why we need Window > Open and Window > Flip command. Also 
a colleague asked me how to demonstrate the behaviour of the sine function sin(x + t) in de- 
pendence of time parameter t. 

Well, the head line is promising too much. I will not describe any ingenious technique in DERIVE and 
a proper animation to be the substitute for AcroSpin or similar program. But for the task described it 
was quite successful. 

The idea is to create the sequence of graphs and by switching simulating an animation. We exploit the 
possibility that DERIVE can open many windows on the same position (Window > Open) and switch 
from one window to the next one (Window > Flip). 

Suppose we like to show how sine waves sin(x + t) “travel” in time. First we prepare a suitable se- 
quence of functions - each one as an own expression: 

SIN (x) 

SIN (x + n/6) 

SIN (x + n/3) 

SIN (x + n/2) 

SIN (x + 1 In/ 6 ) 

Let us plot the last graph: 

Plot > Under 10 > Plot 
Open a new window 

Window > Open > 2D Plot 

Go to the Algebra Window, mark the 
before last expression and plot it 

Algebra 

Plot 

Plot 

We open a new window on the same place and repeat the procedure until we plot the function sin(x). 
Now we have prepared all the pictures and by rapid switching we animate them. Using Window > Flip 
would be too slow, so we use F2. If our computer is sufficiently fast then we will get the impression of 
a travelling wave. 

ACDC - Amusing Corner of the DERIVER' s Curiosity (founded by Alfonso Poblacion Saez) 



It really works!! 





